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CHAPTER  I 

INTRODUCTION 

The  objective  of  this  study  is  to  determine  the 
extent  to  which  infrared  techniques  measure  the  actual 
amplitude  distribution  of  induced  surface  currents  on 
flat  plates  with  given  surface  preparations.  Measure¬ 
ments  of  currents  induced  on  complex  shapes  by  time 
varying  electromagnetic  fields  is,  in  general,  a  time- 
consuming  and  costly  job.  The  reason  is  that  many 
point-by-point  current  measurements  must  be  taken  in 
order  to  characterise  the  induced  current  distribution 
on  a  given  target.  Yet,  knowledge  of  such  current 
distributions  is  essential  to  many  research  programs 
such  as  electromagnetic  pulse  interactions  with 
aircraft  systems,  radar  cross-section  optimization, 
etc . 

The  importance  of  knowledge  of  surface  currents 
in  electromagnetic  problems  can  be  illustrated  with  a 
simplified  approach  to  finding  the  electric  (E)  and 
magnetic  (H)  fields  starting  from  Maxwell's  Equations. 


Ampere's  Law 

VXH-eE  +  J 

(1) 

Faraday ' s  Law 

V  X  E  =  -uH 

(2) 

Gauss'  Law  for 

Electric  Fields  v*E  =  — 

(3) 

e 


Gauss'  Law  —  Magnetic  Fields  =  0  (4) 

The  relationship  between  current  density  (J)  and  charge 


density  (p)  is  found  in  the  equation  of  continuity: 

V*^  “  “it"  (5) 

From  equation  (4)  ,  the  vector  relation  of  the 
divergence  of  a  curl  of  vector  is  equal  to  zero 
(^•vxA  -  0),  and  the  constituent  relation,  B  =  uH,  a 
vector  potential  A  can  be  defined  as  follows: 

ull  =  VXA  (6) 

Now,  substituting  equation  (6)  into  (2)  we  obtain 

VXE  =  -  ~  (VXA)  (7) 

With  the  spatial  and  time  coordinates  treated  independ¬ 
ently,  equation  (7)  can  be  rearranged  as 

3x(g+X)  «  0  (8) 

From  the  vector  identity,  curl  of  a  gradient  of  a 
scalar  is  zero,  from  equation  (8)  one  defines  a  scalar 
potential  V  as 

E+A  -  -VV  (9) 


Rearranging , 

E  =  -VV- A  (10) 

Thus,  the  electrrc  and  magnetic  fields  are  expressed 
in  terms  oi  potentials.  We  can  now  derive  differential 
equations  for  the  potential  function  in  terms  oi  the 
sources  J  and  p . 

Substituting  (6)  and  (10)  into  equation  (1),  we  have 

(VXVXA)  =  —  t  VV— c  A+J  (ID 


Applying  the  vector  identity  ^X^xX  -  to 

equation  (11)  yields 

V2£  -y  t A  «=  -uJ+ut_VV  +  VV  *A  (12) 

Substituting  (10)  into  (3) , 

V2V  +  V  *A  =  -  -J-  (13) 

Now,  equations  (11)  and  (13)  can  be  decoupled  using 
the  Lorentz  gauge  condition. 

-4  -4  • 

V ♦ A  -  -ucV  (14) 


Wave  equations  in  terms  of  a  single  potential  and  source 
terms  follows  the  substitution  of  (14)  into  (11)  and 
(13) 

-4  -4  -4  -4 

V2A  -  ulA  =  -  y J  (15) 


-4  *•  « 

V2V  -  p  tV  =  - 

e 

Finally,  the  solutions  to  the  wa^e  equations  are  of 
the  form 


A  ( r ,  t) 


u 

4*" 


J(r,  t- 
— 


dv  ’ 


(16) 


V  (r ,  t) 


dv' 


(17) 


-4  4  | 

where  r  locates  the  field  point  of  interest  r  locates 
a  source  element,  R  is  the  distance  between  source 
element  and  field  point,  t  -  —  is  the  retarded  time 
correcting  for  the  finite  velocity,  v,  of  the  electro¬ 
magnetic  wave,  and  the  integration  is  performed  over 
the  volume  containing  the  source  lements. 


Figure  1. 

Defining  the  position  vectors. 

Therefore,  knowing  the  magnitude,  direction  and  phase 
of  the  vector  quantity  J  at  all  source  points,  we  can, 
in  principle,  calculate  the  vector  potential  A.  Know¬ 
ing  the  vector  potential,  the  scalar  potential  can  be 
found  through  the  Lorentz  condition,  equation  (14) . 

With  the  vector  and  scalar  potentials,  we  can  now  solve 
for  the  electric  and  magnetic  fields  at  an  arbitrary 
point  in  space.  In  this  study,  only  the  magnitude  of 
microwave  induced  current  density  on  flat  plates  will 
be  investigated.  At  the  microwave  frequencies  used 
in  this  work  the  current  density  is  confined  near  the 
surface  for  highly  conductive  objects  and  is  generally 


called  a  surface  current. 
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This  report  represents  a  necessary  step  toward 

refining  a  technique  which  shows  much  promise  in 

characterizing  induced  surface  current  patterns  with  a 

single  measurement.  The  technique  is  based  on  the  fact 

that  when  current  flow's  in  a  conductor  of  finite 
2 

resistivity,  I  R  heating  occurs.  When  the  surface 
temperature  rises  to  roughly  0.1°  K  or  higher  above 
ambient,  the  induced  temperature  distribution  over 
the  surface  (v’hich  corresponds  to  induced  currents  in 
that  surface),  can  be  detected  by  a  thermoqraohic  system 
via  emitted  thermal  radiation. 

Thermography,  as  applied  to  this  study,  is  the 
heat  detection  process  displayed  in  both  a  color  code 
representation  and  a  numeric  array  of  intensity  values. 
The  science  of  infrared  detection  has  been  advanced 
over  the  past  two  decades  by  the  AGA  Corporation  in 
Sweden.  The  basis  of  operation  is  a  scanning  scheme 
of  infrared  prisms  and  a  photovoltaic  cell  from  which 
voltage  inputs  are  processed.  Applications  of  thermo’- 
graphy  have  spanned  the  areas  of  medicine,  building 
construction,  transmission  line  problems,  blast 
furnace  operations,  etc. 
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At  the  Naval  Postgraduate  School,  Burton  and 

associates  advanced  the  horizon  of  applications  to 

include  the  detection  of  surface  currents  produced  by 

incident  electromagnetic  plane  waves  on  objects,  and 

by  observing  heating  patterns  on  resistive  paper 

produced  by  antennas  (1) .  The  incident  electromagnetic 

wave,  interacting  with  objects  of  appropriate  size, 

produces  standing  wave  patterns,  detectable  in  the 

2 

infrared  region  due  to  the  I  r  heating  effects.  In 
the  laboratory  environment,  the  qualitative  comparison 
of  current  distributions  present  on  antennas  of  various 
configurations  were  also  made  (2) .  Further  results 
demonstrated  a  good  qualitative  comparison  of  induced 
currents  on  flat  plates  and  cylinders  with  values 
measured  at  Harvard.  The  results  appeared  in  a  color 
photograph  representation  such  that  ten  colors  are 
displayed  corresponding  to  a  temperature  range  whereby 
each  color  represented  roughly  one-tenth  of  that  range. 
The  results  were  encouraging  and  would  obviously  be  of 
benefit  if  applied  to  complex  structures  such  as  air¬ 
craft  that  could  be  scaled  down  for  laboratory  examina¬ 
tion.  This  work  indicated  that  a  necessary  step  in 
the  extension  of  the  IR  results  to  various  applications 
was  the  quantifying  of  the  obtained  data.  Once  this 
data  was  obtained  in  this  more  workable  format,  an 
examination  to  determine,  as  accurately  as  possible, 
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the  correspondence  between  actual  surface  currents  and 
those  determined  with  IR  techniques  was  needed.  For 
simple  geometries,  currents  have  been  calculated 
analytically  (22,26)  or  utilizing  numerical  techniques 
(23),  and  verified  experimentally.  For  the  comparison 
study,  square  and  circular  flat  plates  were  logical 
candidates  owing  to  their  high  symmetry. 

The  theoretical  determination  of  surface 
currents  induced  on  the  perfect  electric  conducting 
flat  plates  with  normally  incident  plane  electro¬ 
magnetic  wave  has  been  developed  by  previous  investiga¬ 
tors.  Results  for  circular  flat  plates  of  electrical 
size  of  largest  dimension  .38,  .50,  1.0,  and  1.3 
wavelengths  are  presented.  The  analytic  solution  for 
the  electrically  small  circular  flat  plate  (disk)  will 
be  outlined.  The  perfectly  conducting,  vanishing  thin 
disk  is  computed  through  the  first  six  terms  for  the 
.15  wavelength  disk  and  the  results  are  presented  and 
compared  to  the  solution  obtained  from  an  E-field 
solution  for  a  conducting  body  of  revolution  developed 
by  Joseph  R.  Mautz  and  Roger  F.  Harrington  of  Syracuse 
University  (33).  The  3.0  wavelength  soli  ion  from  the 
Body  of  Revolution  (BOR)  Code  is  also  compared  with  an 
analytic  result  obtained  by  Andrejewski  (  2  ) .  Disks 
of  electrical  size  .38,  .50,  1.0,  and  1.3  are  then 
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solved  using  this  computer  code.  The  flat  plate 
solution  using  a  body  of  translation  electromagnetic 
computer  code  developed  for  the  Rome  Air  Development 
Center  is  presented  and  compared  with  results  from  the 
literature. 

An  infrared/microwave  facility  was  designed, 
constructed,  and  equipped  in  a  converted  storage  room 
at  the  United  States  Air  Force  Academy.  The  available 
microwave  apparatus  dictated  the  size  of  anechoic 
chamber  and  the  frequencies  of  2.45  GHz  and  937  MHz. 
Investigations  of  chamber  characteristics  using 
various  reflector  arrangements  with  a  quarter  wave 
monopole  on  a  ground  plane  were  carried  out  using  a 
ripple-tank  two  dimensional  modeling  apparatus,  a 
crossed-dipole  power  probe,  and  a  time-varying  magnetic 
field  detector.  Corner  reflectors  and  cylindrical 
parabolic  reflectors  were  experimentally  positioned  to 

direct  maximum  power  while  still  achieving  a  "plane" 
wave  over  a  sixteen  centimeter  target.  The  nucleus  of 
infrared  detection  equipment  was  a  AGA  Thermovision 
680  System.  An  interface  for  digitization  was 
developed  such  that  correction  to  the  received  Thermo¬ 
vision  data  could  be  applied  on  a  point  by  point  basis. 
An  infrared  measurement  technique  including  a  seven 
frame  average  of  the  illuminated  matrices  to  reduce 
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random  error  influences  on  the  final  result  was 
developed.  The  averaged  matrix  in  ISU  units  was  then 
corrected  for  camera  calibration  resulting  in  a 
determination  of  temperature  increase  due  to  the 
incident  microwave.  To  obtain  surface  currents  from 
temperature,  a  functional  relationship  was  empirically 
determined  by  applying  various  levels  of  DC  current 
through  the  coating  and  recording  the  corresponding 
temperature  rise. 

The  verification  of  surface  currents  on  metal 
plates  illuminated  with  an  optimized  plane  electro¬ 
magnetic  wave  required  the  design  and  fabrication  of  a 
loop  detector.  The  induced  voltage  is  related  to  the 
magnitude  of  the  time  varying  magnetic  field  through 
Faraday's  Law  and  hence  to  surface  current.  Considera¬ 
tion  of  shielding  from  stray  RF  coupling  as  well  as 
electric  field  coupling  were  found  to  be  critical  to 
obtain  accurate  measured  values. 

Finally,  current  values  obtained  from  theory, 
probe  measurement,  and  IR  detection  are  compared.  The 
comparison  procedure  developed  should  be  reapplied  for 
future  coating  schemes  using  higher  resolution  infrared 
detection  systems.  Data  gained  through  correlation 
studies  on  simple  shapes  will  then  permit  an 
extension  to  more  complex  shapes. 


Supplemental  studies  in  wave  interaction 
modeling  and  emissivity  dependence  on  angle  were  also 
conducted  and  are  summarized  in  the  Appendices. 
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CHAPTER  II 

THEORETICAL  DETERMINATION  OF  CURRENTS  ON 
FLAT  PLATES 

The  solution  to  the  metallic  plate  or  even  the 
perfectly  conducting  plate  is  not  trivial  as  emphasized 
in  a  paper  given  by  Professor  T.  B.  A.  Senior  in  May 
1982  (45)  in  which  recent  approaches  to  solve  the 
integro-differential  equations  for  arbitrary  flat 
plates  given  in  a  1977  IEEE  publication  were  found  to 
be  in  error  as  the  plate  shape  varied  from  the  elliptic 
case.  For  perfectly  conducting  disks  of  largest 
largest  electrical  dimension  .38,  .50,  1.0,  and  1.3 
wavelengths,  a  Body  of  Revolution  electromagnetics  code 
was  used  for  their  respective  solution  at  normal 
incidence.  Since  the  disk  has  been  solved  analytically, 
the  numerical  approach  of  the  computer  code  was  com¬ 
pared,  at  both  extremes  of  the  experimental  dimensions, 
to  an  analytical  treatment.  Analytic  solutions 
developed  by  C.  J.  Bouwkamp  (7)  and  W.  Andrejewski  (2) 
were  used  for  the  .15  and  3.0  wavelength  circular  disks, 
respectively.  The  code  results  showed  good  agreement 
with  the  analytic  solutions.  For  the  perfectly 
conducting  square  plates,  solutions  from  the  literature 


12 


are  presented  in  the  range  of  the  experimental 
sizes.  An  additional  Body  of  Translation  electro¬ 
magnetics  computer  code  solution  is  found  in  Appendix 
C  for  the  plate  sizes  used  in  this  report. 


Circular  Flat  plate  (Disk) 


The  general  formulation  for  this  problem 

begins  with  the  requirement  that  the  electromagnetic 

field  must  obey  Maxwell's  equations.  At  the 

surface,  the  boundary  conditions  must  be  satisfied. 

In  the  perfectly  conducting  case,  the  tangential 

compo.  i  of  the  electric  field  and  the  normal 

compone*  of  the  magnetic  field  must  vanish.  For 

flat  objects  the  special  situation  at  the  edges  must 

also  be  coi.  idered.  If  one  observes  that  current  or 

charge  ca;.  .eside  at  an  edge  and  that  the  electric 

field  varies  as  1/r,  one  at  first  would  allow,  on  a 

theoretical  basis,  for  the  electric  fields  and  current 

to  approach  infinity  at  the  edge.  A  generally 

accepted  conservation  of  energy  argument  dealing  with 

the  ord  r  of  this  mathematical  singularity  is  stated 

by  Josef  Meixner  as  follows  (34): 

The  edge  condition  states  that  the 
electromagnetic  energy  density  must 
be  integrable  over  any  finite  domain 
even  if  this  domain  contains  sing¬ 
ularities  of  the  electromagnetic 
field.  In  other  words,  the  electro¬ 
magnetic  energy  in  any  finite  domain 
must  3e  finite.  In  the  case  of  a 


perfectly  conducting  surface 
with  an  edge,  one  concludes  from 
this  condition  that  near  the  edge 
the  singular  components  of  the 
electric  and  magnetic  field 
vectors  are  of  the  order  p"1'2, 
where  p  is  the  distance  from  the 
edge,  while  the  components  of  the 
electric  and  magnetic  field 
strengths  parallel  to  the  edge 
are  always  finite. 

The  following  discussion  is  taken  from 
Bouwkamp  (7)  to  illustrate  the  solution  process  for 
the  simplest  case,  that  of  an  electrically  small, 
perfectly  conducting  circular  disk  illuminated  by  a 
plane  wave  at  normal  incidence.  The  disk  is  taken  to 
be  of  unit  radius  and  centered  on  the  origin  of 
rectangular  coordinates  X,  Y,  Z  so  that  the  rims  of 
the  disk  is  defined  by  Z  =  0,  X2+Y2  =  1  (see  Figure  2) 


Figure  2.  Coordinate  system  for  the  disk. 


The  incident  wave  is  plane  polarized  with  the  electric 
field  along  the  X  axis,  the  magnetic  field  aligned 
with  the  Y  axis,  and  thus  the  wave  is  moving  in  the 
positive  Z  direction.  This  solution  is  restricted  to 
ka<l  where  k  is  the  wave  number  and  a  is  the  radius  of 
the  disk.  The  general  expressions  relating  the  field 
vectors  with  the  scalar  and  vector  potentials  (equa¬ 
tions  6  and  10)  and  the  general  expression  for  the 
vector  potential  in  terms  of  current  density  serve  as  a 
starting  point.  The  scattered  field  can  then  be 
written  in  rectangular  coordinates  in  terms  of  the 
vector  potential  as 
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The  boundary  conditions  for  a  perfectly  con¬ 
ducting,  and  thus  perfectly  reflecting,  disk  are  that 
the  total  tangential  electric  field  and  the  total 
normal  magnetic  field  must  vanish  at  the  surface.  The 
total  field  is  simply  the  incident  plus  the  scattered 
field.  Thus,  at  the  surface  (X2+Y2-l) 

Ex  =  “  Ex  <24> 

Ey  =  0  (25) 

Hz  =  0.  (26) 

Now  from  the  scattered  field  expressions  and  the 
stated  boundary  conditions,  the  following  simultaneous 
system  of  equations  can  be  written  as  follows: 


3 Ay  _  3  A 


9y  3x 

3  2  Ax  3  2  Ax 
3x2  3y2 


+ 


k2Ax 


+  k2a. 

3x2  3y2 


ik 


0 


with 


1 

Ay  ~  c 


Jy 


ikr 

e 


r 

ikr 

e 


r 


ds 


ds 


(27) 

(28) 

(29) 

(30) 

(31) 


These  are  the  integro-dif ferential  equations  that  must 
be  solved  to  obtain  the  rectangular  components  of  the 


L6 


current  density.  Now,  Bouwkamp  reduces  this 
simultaneous  system  to  that  of  an  ordinary  system  by 
using  symmetry  and  normal  incidence  argument  to 
eliminate  the  functions  and  A^.  What  is  left  is 
set  of  integro-dif ferential  equations  in  only  Jx. 

The  conditions  on  the  edge  of  the  disk  must  be 
considered.  From  a  combination  of  physical  and 
mathematical  arguments  stated  earlier,  order  functions 
are  obtained.  These  state,  that  near  the  rim,  as  the 


distance  to  the  rim  goes  to  zero,  the  current 
component  on  the  surface  of  the  disk  that  is  normal  to 
the  edge  will  decrease  with  the  square  of  the  distance 
and  the  current  component  on  the  surface  that  is 
tangential  to  the  edge  will  increase  as  the  reciprocal 
of  the  square  root  of  the  distance. 

To  take  advantage  of  the  symmetry  in  this 
situation,  the  incident  wave  is  considered  in  terms 
of  cylindrical  coordinates.  Functions  depending  on 
the  radius  (p)  are  then  introduced  along  with 


(32) 
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where  J  and  J  are  related  to  each  other  by 
x  y 


j  =  _  -  9fx  =  k  (y  i£x  _x  9 
Jy  2  2  ^  ax  X  3y 


(33) 


The  functions  A  and  B  must  also  satisfy  the  conditions 
that  on  the  rim  their  sum  is  zero  and  at  the  center  of 
the  disk  the  value  of  B  is  zero. 

Now,  solutions  for  equations  (31)  and  (32) 
eventually  reduce  to  being  able  to  solve  the  following 
integral  equation: 


p  'dp  ' 


A ( p  1 )  +  B  (p  ')  cos2<t> ' 
it2  y i-p ^ 


-  -  P  q,  (kp)  +  p  J?  (kp)  cos2* 

for  3  unknowns  p,  A,  B 
over  oip-1  and  o~4>-2tt 


(.34) 


where  r2  =  p2  -  2pp*  cos(<f-if>')  +  p'2 

A ( p ' 1  and  B(p')  are  uniformly  bounded  on 
the  internal  and  JQ  and  J2  are  Bessel  functions 
of  argument  kp  and  of  order  0  and  2  respectively. 
From  a  power  series  solution  of  equation  (34) 
are  found  series  solutions  of  A ( p )  and  B(p)  from 

equations  (31)  and  (32).  For  the  first  six  terms, 

6 

A(P)  =  E  An(ik)n 
n=l 

=  j-  (-4+3p2).ik  +n0  + 

+  i  (56-40p2+5p4)  (ik)3 


+  —  (2-pz)  (ik) 14  + 
9 


+  rr^rr  (-2656  +  24Q8p2  -  448p4  +  21p6)  (ik) 5 
12600 


+  ~  (-296  +  192p2  -  15p4)  (ik)6  +  O  (A?) 


and  B(p)  =  ^  Bn  (ik) 
n=l 


=  3  P  2  (  ik)  +  Ott  + 

+  p2  ( -8  +  p 2 )  (ik) 3  + 

-  y-  p2  (ik)4  + 

+  75-20  °2  (  200-68p  2  +  3p  4 )  (ik)5  + 

+  ™  p2  ( 1 34-15p 2 )  (ik)6  +  0  (B;) 


Simplifying, 

A  ( o )  =  k 


12600 

2jnk5 

675 


B  (  p )  =  k 


k  {  3  [-« + 3  <£>  ]  + 

J  56-40 (|) 2  +  5(|) “ ] 

[  2  -  <!>  J 

■  £  -2656  +  2408  (|)  2  -  488  (|)  -  21(|) 

[  -  296  +  192  (|)  2  -  15  (|)4  ]  +  Oj  ( 

<i»2  -¥  ‘S>2  [  -«  -  <!>2  J 


4  nk  2 
9 


(&)  + 


'a 


(P)2  f  200-68  (|)2  +  3(|) 


kipj 
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iiJd  (P)2 

675 


134  -  15  (£) 


]+oj 


(38) 


Now,  Bouwkamp  assumes  in  the  initial  development  that 

e  y  u 

=  1.  However,  g-  =/  ^°  =  377  and  therefore, 

1  o 

Jx,  Jy  must  include  -jyy  factor. 


To  find  the  magnitude  of  the  current  components  we 
will  simplify  by  separating  the  real  and  imaginary 
parts 

i  i  -  /J-.<  jJ 


(39) 


J 


ck 


377tt  y/  1-(P) 
[A(p)  +B(p)  cos2<*] 


£[a(p) +B(p)  cos2<J>  J 


(40) 


J.  = 


k5 

675' 


ck  (47.  v  r  o_/P\2*|_  2tt 

X  K  L  2  j  ~6 

3777.2/ l-(f)  2  C 

cl 

[-296+192  (|)2-15  (s)^+j0  [-4+3  (J)2]  - 


90 


[56-40  (§)  +  5  (§)  j  +  [-2656  +  2408  (|) 


D  4  o  6 

-  448  (|)  +  21  (£) 


115 


kJ  ,P 


2”k5  ,£. 
675  la; 


[134  -  15  (§>*  ]J  cos2^  +  j  COS2*  ^|  [  (|)  2  ] 

111  (Pi2r-a  +  f  P  x  2l  .  k4  ,p  2  [200-68  (J)  2 
”  30  L  8  +  (a>  -1  +  2520 

+  3  (|)4  jj  (41) 


Jx  ’  - —f= f 

377,2^  /  1-(|) 


where 


[a^  +  ja2)  +  (b1cos24i  +  ]b2 

cos2$)]  (42) 


=  real  part  of  A(p) 
a2  =  imaginary  part  of  A(p) 
b-^  =  real  part  of  B(p) 
b2  =  imaginary  part  of  B  ( p ) 

O 

Simplifying  further,  c  =  3.00  x  10  m/sec 
,  _  ,  5.05  x  105  Cr 

lJxl  I  p  2  7  L(ai+b1cos2<{)  +  j  (a2+b2cos2$] 

*  C 

[(a^+b^cos2<()  )  -  j  (a2+b2COs2<|>)3P 
5 

-  5-‘-° [  (a^+b2cos24>)  2  +  (a2+b2cos2<ji)  2J  % 

I  /s  2  f  A 


X  /  1-(|) 


Similarly , 


IV  =  (JyJy*>15 


_ c2ti  sin2i{) 

377ti  2  X  /  1-(J) 


5.05  x  IQ3  sin2( 

x  / 1-(|>2 


(B  (  p )  BCp)*)^ 


(bl2  +  b22)JS 


Now,  let  the  total  current  be  Jt  which  is  given  by  the 
following  expression: 


2 


2 


5 


Consider  the  term 


in  the  denominator 


of  Jx  and  Jy.  As  p  approaches  the  radius  a,  the 
currents  tend  to  infinity.  The  calculations  based 
on  equations  (  44  )  and  {  45  )  will  be  restricted  to 

p_^.95a. 

It  was  advantageous  to  transform  from 
cylindrical  coordinates  to  rectangular  coordinates 
using  the  following  substitutions: 


sin2$ 

cos2<l) 


2  tan  (j> 

1  +  tan2<{> 

1  -  tan2<{> 
1  +  tan2<|> 


where  tan  <f  =  x 


It  follows  that 


sin24>  = 


cos24>  = 


2xy 


x2+y2 


x2_ y2 


x2+y2 


The  latter  equations  were  then  programmed  on  an 
HP9845  minicomputer  and  the  results  for  a  .075 A  radius 
disk  are  given.  With  a  vertically  polarized  incident 
plane  wave,  the  total  current  density  along  the 
horizontal  and  vertical  axis  can  be  computed.  Curve 
fits  for  the  results  follow  a  discussion  of  the 


program  used. 
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The  regression  curve  fits  were  done  using  a 
Hewlett  Packard  software  package  on  the  9845B  mini¬ 
computer.  The  program  takes  a  set  of  points  (X,  Y) 
and  fits  various  types  of  curves  to  them  using  least 
square  methods.  An  analysis  of  variance  table  is  also 
computed  and  this  information  is  printed  in  the  text 
that  follows.  The  input  parameters  (X^,  Y^ )  are 
loaded  into  full-precision  one-dimensional  arrays  of 
N  elements  each  with  being  at  least  3.  The  output 
variables  are  also  full  precision  and  are  defined  as 
follows : 


Regression 

Df  -  regression  degrees  of  freedom 
(degree  for  polynomial  model) 
SS  -  regression  sum  of  squares 
MS  -  regression  mean  squares 
Residual 


Df  -  residual  degrees  of  freedom  (N-2) 
SS  -  residual  sum  of  squares 
MS  -  residual  mean  squares 


Total 

Df  -  total  degrees  of  freedom  (regression 
and  residual) 

SS  -  total  sum  of  squares 
F  -  the  F  statistic  (regression  mean  squares/ 


residual  mean  squares) 


The  F  ratio  is  an  indication  of  how  good  the 


curve  fit  is  and,  in  general,  is  optimized  by  the 
highest  possible  F  ratio.  In  the  general  outputs 
that  follow,  a  blank  space  for  an  F  ratio  implies 
an  exact  fit  such  that  all  points  fall  on  the 
curve.  It  should  also  be  noted  that  several  factors 
enter  the  goodness  of  fit  such  as  the  number  of  data 
points,  the  number  of  coefficients,  etc.  For  data 
sets  in  this  paper,  only  the  linear  and  polynomial 
models  were  found  to  be  appropriate.  For  n  data 
points  ( )  the  method  and  formulae  used  by  the 
program  are  as  follows  (24): 

For  the  linear  model  y  =  a+bx: 


ZXiYi  -  £x.  ZYj 

b  =  - 2 - 

Zxi2  -  (Ix^2 

n 

Total  sum  of  squares  =  zyi2  -  ( y i ) 2 
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Regression  sum  of  squares  -  [Zx.-y^  -  Exi  ly.^]2 

_ n _ 

£xi2  -  (Ix^2 
n 


Residual  Sum  of  Squares  =  Total  SS  -  REG  SS 

Regression  Mean  Squares  =  REG  SS 

Residual  Mean  Squares  =  RES  SS 

n-2 


F  Ratio 
For  the 

y  = 


RES  MS 
RES  MS 


polynomial  model 


m  ,  m-1 

ax  +  a  1  x 
m  m-1 


+ 


+a i x  +  a  o  : 


m  =  degree  of  the  desired  polynomial 

(x^,  f  (xi) )  =  set  of  tabulated  points  ((Xj,y.j.)  in  the 
other  models 

gj (xi )  =  xiJ  for  i  =  0,  1,  ....  n 


akj  =  ,E  8k(xi)gj(xi)  for  k  =  °’ 

j  =  0, 

A  system  of  simultaneous  equations  is 
in  terms  of  a's  as  follows 


a  a  +  a  a  +  . . .  +  a0,a(  =  Z 
mom  m- 1 o  m-i  u  u  u 

}  -  o 


a  a  +  a  a  +  ...  +  a0ia0  =  £ 


m  i  m  m- i i  m- i 


i=o 


1 ,  2,  .  .  .  ,  m 
1 ,  2,  .  .  .  ,  m 
then  written 

f(xjL)  g  0  (  X  .  ) 

f(xt)  gi(x.) 
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a  a  +  a  a  + 
m 2  m  m—  1 2  m—  i 


a  a  +  a  a  +  , 
mm  m  m-im  m-i 


+  a0  2a0  =  £  f (x±  )  g2  (xi  ) 

i=  o 


•  +  °.ma"  "  E  f(xi>  gm  <xl> 

1=  o 


which  can  be  solved  to  find  the  coefficients 


a  o  ,  a i ,  a2 i ... i a 


y  = 


yi  *  Vi"1  +  am-.  xi'"'1  +  •••  +  a-xi  +  a« 


Total  Sum  of  Squares  =  £  (y^-  y) 


Regression  Sum  of  Squares  =  Z(yi-  y)2 
Residual  Sum  of  Squares  =  Total  SS  -  REG  SS 


Regression  Mean  Squares  =  REG  SS 


Residual  Mean  Squares  =  RES  SS 

n-l-m 
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DATA 


Po  i  nt 

#1 : 

X  ™  0 

Y=5 . 45 

Point 

#2: 

X=1 

Y-5.41 

Point 

#3: 

X  =  2 

Y=5. 46 

Poi  nt 

*4: 

X  =  3 

Y*5. 34 

Poi  nt 

#5: 

X  =  4 

Y  =  5 . 48 

Point 

#6: 

X*5 

Y=5 . 22 

Po  i  nt 

#7: 

X  =  6 

Y  =  5 . 5 

Po  i  nt 

#8: 

X  =  7 

Y  =  5 . 5 

Po  i  nt 

X  =  8 

Y=5 . 54 

Poi  nt 

#10: 

X  =  9 

Y  =  4 . 83 

Po  i  nt 

#11: 

X=10 

Y=5. 59 

Po  i  nt 

#12: 

X*1 1 

Y*4. 56 

Point 

#13: 

X=  1 2 

Y®5. 67 

Poi  nt 

#14: 

X=  1 3 

Y=4. 21 

Poi  nt 

#15: 

X=  1 4 

Y=5. 81 

Poi  nt 

#16: 

X=  1 5 

Y*3. 77 

Po  i  nt 

#17: 

X=  1 6 

Y=6. 04 

Po  i  nt 

#18: 

X=  1 7 

Y=3. 21 

Po  i  nt 

#19: 

X=  1 8 

Y=6. 47 

Poi  nt 

#20: 

X=  1 9 

Y=2. 44 

Poi  nt 

#21 : 

X  =  20 

Y=7. 46 

Poi  nt 

#22: 

X  =  2 1 

Y= 1 . 16 

Point 

#23: 

22 

Y= 1 1 . 3 

Figure  3:  Data  from  an  analytic  solution 

(Bouwkamp)  for  the  .15  wavelength 
circular  disk. 


—i — i  i  i  i  !  i  i  i  i 

ojnvintONfflms-w 


CO  Tj-  ID  03  N  CO  cn 


Distance  From  Disk  Center  (m) 


Plot  of  Jp  and  J<p  current  components  at 
<P  =  0°  and  <j>  =  90°  for  the  .15  wavelength 
circular  disk  (analytic  solution) . 


20  J 
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DATA 


Poi  nt 

#1 : 

X  =  0 

Y  =  5 

.45 

Poi  nt 

#2: 

X=  1 

Y-5 

.41 

Poi  nt 

#3: 

X  =  2 

Y  =  5 

.34 

Poi  nt 

#4: 

X*3 

Y«5 

.22 

Poi  nt 

#5: 

X  =  4 

Y  =  5 

.05 

Poi  nt 

#6: 

X  =  5 

Y  =  4 

.  83 

Poi  nt 

#7: 

X  =  6 

Y-4 

.56 

Poi  nt 

#8: 

X  =  7 

Y  =  4 

.21 

Poi  nt 

#9: 

X  =  8 

Y  =  3 

.77 

Po  i  nt 

#10: 

X  =  9 

Y*3 

.21 

Poi  nt 

#u: 

X=10 

Y  = 

2.44 

Po  i  nt 

#12: 

X=1 1 

Y  = 

1 .  16 

POLYNOMIAL  MODEl:  Y  =  A<t1>*X''M+A(M-l  >*X"<M-1  )  +  ...+A(l>*X  +  A<0.> 
Coef f i c i entsl 

A <@>=5,454547487 
07842408 
A <2 >=.03991346 
A (3 >=-. 019131686 
R<4>=. 0024157556 
A<5)“-. 0001 1274502 


Source 

Df 

SS 

MS 

Regress i on 

5 

20. 1 13 

4,023 

12702. 7636366 

M4D. 3D 

Res i dual 

6 

.  002 

.  000 

T ot  al 

1 1 

26. 115 

Figure  5:  Polynomial  fit  for  the  Jp  current 

component  at  $  =  0°  for  the  .15 
wavelength  circular  disk  (analytic 
solution) . 


Plot  of  the  Jp  current  component  at 
<t>  =  0°  for  the  .15  wavelength  circular 
disk  (analytic  solution) . 
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DATA 

Point  #  1 : 
Point  #2: 
Point  #3: 
Point  #4: 
Point  #5: 
Point  #6: 
Point  #7: 
Point  *8: 
Point  #9: 
Point  # 10 : 
Point  #11: 


X*0 

Y*5. 45 

X=  1 

Y-5.46 

X  =  2 

Y=5. 48 

X  =  3 

Y=5.5 

Xa4 

Y-5.54 

X  =  5 

Y*5. 59 

X-6 

Y«5. 67 

X»7 

Y*5. 81 

X<*8 

Y=6. 04 

X-9 

Yb6 . 47 

X^10 

Y*7. 46 

POLYNOMIAL  MODEL:  Y=fl < M > *X-M+A< M- 1  ) ♦X'" C M- 1 )  +  .  .  .  +  A <  1) *X+A < 0 > 
Coef f i c i ents : 

A<0>=5. 445332077 
A<1>=. 87520403 
A<2)=-. 06801385 
A<3>=. 02502269 
A<4>=-. 003571543 
A<5>=. 00018750039 


Source 


Df 


SS 


MS 


F 


Regression  5 
Residual  5 
Total  10 


3.779 

.002 

3.780 


.756  2194.413 

.000 


Figure  7:  Polynomial  fit  for  the  J 4>  current 

component  at  $  =  90°  for  the  .15 
wavelength  circular  disk  (analytic 
solution) . 


Plot  for  the  J$  current  component  at 
<j)  =  90°  for  the  .15  wavelength 
circular  disk  (analytic  solution) . 
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Now,  Bouwkamp's  solution  is  limited  to  ka<l 

2  T1  • 

where  k  =  and  a  is  the  radius  of  the  disk.  To 
obtain  theoretical  solutions  on  the  disks  of 
experimental  interest,  a  Body  of  Revolution  electro¬ 
magnetics  computer  code  developed  by  J.  R.  Mautz  and 
R.  F.  Harrington  (33)  was  used  with  the  input  data 
appearing  in  Appendix  B.  In  the  figures  that 
immediately  follow  for  the  .15  wavelength  circular 
plate,  the  correlation  with  the  solution  of 
Bouwkamp  is  good. 

In  all  data  plots  for  the  J 4>  current  component 
it  should  be  noted  that  the  point  nearest  to  the  rim 
departs  significantly  from  the  neighboring  interior 
point.  The  polynomial  modeling  was  unable  to  follow 
this  departure  from  a  relati'.ely  smooth  curve.  The 
last  point  was  therefore  deleted  from  the  polynomial 
fit  to  the  J<f>  component  throughout  this  text.  This 
edge  point  deletions  should  be  recalled  in  Chapter  VII 
that  there  exists  a  very  high  edge  current  associated 
with  the  maximum  value  of  J$,  though  not  depicted. 
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DATA 


Po  i  nt 

tti : 

X*.001 

Y-.8792 

Po  i  nt 

#2: 

X-.002 

Y*. 8768 

Po  i  nt 

#3: 

X* . 003 

Y*. 872 

Po  i  nt 

*4: 

X* . 004 

Y-.8637 

Po  i  nt 

#5: 

X*. 065 

Y*. 8731 

Po  i  nt 

#6: 

X* . 006 

Y-.8464 

Po  i  nt 

#7: 

X-. 007 

Y-.8741 

Po  i  nt 

#8: 

X*. 008 

Y=. 8224 

Po  i  nt 

#9: 

X-. 009 

Y-.8776 

Po  i  nt 

tie: 

X«.  01 

Y-.7911 

Po  i  nt 

#n: 

X».011 

Y* . 8821 

Po  i  nt 

#12: 

X*. 012 

Y*. 7514 

Poi  nt 

#13: 

X*. 013 

Y*. 8923 

Po  i  nt 

#14: 

X*. 014 

Y* . 7021 

Po  i  nt 

#15: 

X*.  0 15 

Y=. 9091 

Po  i  nt 

#16: 

X«.  016 

Y-.641 

Po  i  nt 

#17: 

X*. 017 

Y=. 9409 

Po  i  nt 

#18: 

X3 .018 

Y«. 5645 

Po  i  nt 

#19: 

X*. 019 

Y= 1 . 007 

Po  i  nt 

#20: 
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Figure  10:  Plot  of  Jp  and  J$  current  components  at 

<f>  =  0°  and  $  =  90°  for  the  .15  wavelength 
circular  disk. 
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Figure  11:  Polynomial  fit  for  the  Jp  current 
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wavelength  circular  disk. 
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Figure  13:  Polynomial  fit  for  the  J (J>  current 

component  at  <p  =  90°  for  the  .15 
wavelength  circular  disk. 
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circular  disk. 


To  confirm  that  the  code  is  reasonable  at  the 
longer  wavelengths,  a  3.0  wavelength  disk  was  run  and 
compared  to  the  analytic  solution  of  Andre jewski  (2) 
which  immediately  follows  the  plot  from  the  code  (Figure 
17).  The  correlation  is  again  quite  good.  After  the 
summary  graphs  (Figure  19)  are  found,  the  polynomial 
fits  for  the  maximum  Jp  and  J<j>  current  components  for 
the  experimental  circular  disks  of  diameters  .38, 

.50,  1.0,  and  1.3  wavelengths  are  given. 


Po  i  nt 

#i: 

X= . 0063 

Y=3. 441 

Po  i  nt 

#2: 

X  =  . 0122 

Y=3. 285 

Po  i  nt 

#3: 

X= . 0 1 85 

Y=3. 106 

Po  i  nt 

#4: 

X= . 0245 

Y=2. 66? 

Po  i  nt 

#5: 

X= . 0307 

Y=2. 906 

Po  i  nt 

#6: 

X* . 0367 

Y=l. 953 

Po  i  nt 

#7: 

X=. 0429 

Y=2. 62? 

Po  i  nt 

#8: 

X=. 049 

Y=1 . 29 

Po  i  nt 

#9: 

X=. 055 1 

Y=2. 358 

Po  i  nt 

#10: 

X= . 061 2 

Y=. 8748 

Po  i  nt 

#11 : 

X=. 0673 

Y=2. 125 

Point 

#12: 

X=. 0734 

Y=. 8345 

Point 

#13: 

X= . 0795 

Y=l„ 96? 

Po  i  nt 

#14: 

X*. 0857 

Y=  1 .  149 

Po  i  nt 

#15: 

X=. 0917 

Y=  1 . 886 

Po  i  nt 

#16: 

X= . 0979 

Y=1 .706 

Po  i  nt 

#17: 

X=. 1039 

Y-l. 875 

Po  i  nt 

#18: 

X  = .  1 1.02 

Y=2. 343 

Po  i  nt 

#19: 

X= . 1161 

Y=l. 899 

Po  i  nt 

#20: 

X= . 1224 

Y=2. 869 

Po  i  nt 

#2i: 

X= . 1283 

Y= 1 . 925 

Po  i  nt 

#22: 

X= . 1346 

Y=3.  115 

Po  i  nt 

#23: 

X=. 1405 

Y=  1 . 918 

Po  i  nt 

#24: 

X=. 1469 

Y»2. 989 

Point 

#25: 

X= . 1527 

Y  =  1 . 866 

Po  i  nt 

#26: 

X=. 1591 

Y=2. 4585 

Po  1  nt 

#27: 

X=. 1649 

Y=1 . 778 

Po  i  nt 

#28: 

X=. 1714 

Y=1 . 667 

Po  i  nt 

#29: 

X=. 1771 

Y=2. 898 

Figure  15:  Data  from  the  Body  of  Revolution 

Code  for  the  3.0  wavelength 
circular  disk. 
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Figure  16:  Plot  of  Jp  and  J<p  current  components  at 

<}>  =  0 0  and  <j>  =  90°  on  the  3.0  wavelength 
circular  disk. 
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Figure  17:  Theoretical  result  for  3.0 

wavelength  circular  disk  (  2  )  . 


In  this  figure  is  shown  the  illuminated  side  of 
a  disk  of  approximate  diameter  of  3.0  wavelengths.  The 
top  graph  corresponds  to  Jp  at  <t>  =  0°  and  the  left  graph 
corresponds  to  J<j>  at  <J>  =  90°. 


In  Figure  19  are  summary  graphs  for  the  perfectly 
conducting  disk  in  which  C  =  ka.  Note  that  at  approx¬ 
imately  1.2  wavelengths  in  diameter  that  the  Jp  component 
begins  to  depart  from  the  case  of  the  small  disk. 
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Figure  20:  Data  from  the  Body  of  the  Revolution 

Code  for  the  .38  wavelenth  circular 
disk . 


Po  i  nt 

#1 : 

X=. 0051 

Y=3. 327 

Po  i  nt 

#2: 

X=. 0101 

Y=3. 27 

Po  i  nt 

#3: 

X  = .  0 1  5  2 

Y=3. 191 

Po  i  nt 

#4: 

X=. 0203 

Y=3. 084 

Po  i  nt 

#5: 

X*. 0253 

Y  =  2 . 9  5 

Po  i  nt. 

#6: 

X=. 0304 

Y=2. 78 

Po  i  nt 

#?: 

X=. 0355 

Y  =  2 . 5  7  2 

Po  i  nt 

48: 

X=. 0405 

Y=2. 326 

Po  i  nt 

#9: 

X=. 0456 

Y-2. 023 

Po  i  nt 

#10: 

X=. 0507 

Y= 1 . 643 

Point. 

#11 : 

X= . 0557 

Y«l. 166 

POLYNOMIAL  MODEL:  Y=A<N>*X  VM+A<M-1  >*X/'<M-1  )  +  .  .  .+AC1>*X+A(0> 
Coefficients: 

A <0> =3. 402624814 
AC  1 >=-13. 5166622 
A<2>=87. 090219 
Ac3> =-10076, 8948 


S  ourc e 

Df 

ss 

MS 

Regressi on 

3 

5.089 

1 . 696 

13340.7291375 

M4D. 3D 

Residual 

7 

.  001 

.000 

T  o  t  a  1 

10 

5.090 

Figure  21: 

Polynomial 

fit  for  the 

Jp 

current 

component 

at  <t>  =  0°  on 

the 

.38  wave 

length  circular  disk. 

Current  Den 


Figure  22: 


Plot  of  the  Jp  current  component  at 
=  0 0  for  the  .38  wavelength  circular 
disk  . 
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Figure  23:  Polynomial  fit  for  the  current 

component  at  <j>  =  90°  for  the  .38 
wavelength  circular  disk. 
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Figure  26: 


Plot  of  Jp  and  J<t>  current  components  at 
<fi  =  0°  and  $  =  90°  on  the  .50  wavelength 
circular  disk. 
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Figure  27:  Polynomial  fit  for  the  Jp  current 

component  at  <j>  =  0°  for  the  .50 
wavelength  circular  disk. 
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Figure  28:  Plot  of  the  Jp  current  component  at 

<p  =  0°  for  the  .50  wavelength  circular 
disk . 


Plot  for  the  Jif>  current  component  at 
$  =  90°  for  the  .50  wavelength 
circular  disk. 
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Code  for  the  1.0  wavelength 
circular  disk. 
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Figure  32:  Plot  of  Jp  and  J 4>  current  components  at 

<f>  “  0°  and  $  =  90°  for  the  1.0  wavelength 
circular  disk. 


rrent  Densitv 


Figure  34:  Plot  of  the  Jp  current  component  at 

<f>  =  0°  for  the  1.0  wavelength  circular 
disk  . 
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component  at  <(>  =  90°  for  the  1.0 
wavelength  circular  disk. 
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Data  from  the  Body  of  Revolution 
Code  for  the  1.3  wavelength 
circular  disk. 
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Figure  39:  Polynomial  fit  for  the  Jp  current 

component  at  <t>  =  90°  for  the  1.3 
wavelength  circular  disk. 
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<f>  =  0°  for  the  1,3  wavelength  circular 
disk. 
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Figure  42:  Plot  for  the  J<f>  current  component  at 

$  =  90°  for  the  1.3  wavelength 
circular  disk. 
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The  square  flat  plate  is  not  solvable 
analytically.  A  Body  of  Translation  electromagnetic 
computer  code  solution  was  run  at  the  Rome  Air 
Development  Center  and  data  results  were  then  collated 
into  tables  for  total  current  for  the  top  half  of 
the  square  plate  and  contour  plots  were  drawn.  These 
appear  in  Appendix  C  along  with  the  current  components 
at  the  horizontal  bisecting  line  of  a  1.0  wavelength 
square  plate.  This  data  differs  from  that  of  other 
codes  such  as  the  ones  that  appear  in  the  next  three 
figures.  The  square  plates  are  of  electrical  sizes 
.15,  1.0,  and  1.2  wavelengths  and  the  dominant  current 
component  is  vertical  if  we  consider  the  incident  El- 
field  to  also  be  oriented  vertically.  It  should  be 
noted  that  the  plates  are  assumed  to  be  perfectly 
conducting  and  that  the  form  of  the  current  density 
plots  do  not  vary  significantly  from  .15  to  1.2  wave¬ 
lengths.  Only  a  slight  rise  in  current  is  observed 
near  the  midpoint  of  1.0  and  1.2  wavelength  plate 
unlike  the  case  of  the  .15  wavelength  plate.  Also, 
shown  in  Figure  45  is  the  current  on  a  resistive  plate 
depicted  by  a  dashed  line. 
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Figure  45:  Theoretical  result  for  1.2  wavelength 

square  plate  (45). 
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CHAPTER  III 

•  INFRARED/MICROWAVE  LABORATORY 

To  conduct  an  experiment  for  detecting  micro- 
wave  induced  surface  current  through  an  infrared 
approach,  it  was  necessary  to  design  and  construct  a 
facility  where  none  had  previously  existed.  For  a 
"proof  of  principle"  goal,  flexibility  was  highly 
valued  as  unforeseen  problems  were  bound  to  occur. 

Time  and  financial  constraints  dictated  that  available 
resources  be  used  to  the  maximum  extent  possible.  To 
observe  heating  effects  on  relatively  preliminary 
coating  schemes,  the  highest  microwave  levels 
possible  from  available  sources  were  to  be  obtained, 
thus  safety  considerations  were  also  paramount.  The 
basic  experimental  arrangement  is  seen  on  the  follow¬ 
ing  page  (Figure  46). 

The  antenna,  metal  parabolic  and  corner 
reflectors,  and  target  object  were  mounted  on  a  12'  x 
4'  aluminum  ground  plane  which  was  encased  by  absorb¬ 
ing  cones  forming  an  anechoic  chamber.  The  electro¬ 
magnetic  wave  propagates  from  the  antenna/reflector 
configuration  to  the  target.  The  wave  interaction 
with  the  target  creates  surface  currents  which  in  turn 
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produce  joule  heating  detectable  by  the  Thermovision 
camera  located  behind  the  reflectors.  The  camera 
was  maintained  at  normal  incidence  from  the  target 
surface. 

The  basic  setup  is  seen  in  the  following 
photograph  (Figure  47) .  Located  to  the  left  of  the 
chamber  on  a  cart  is  the  apparatus  for  recording 
position  and  induced  voltage  from  a  magnetic  field 
probe  which  was  developed  to  detect  surface  currents. 
The  magnetic  field  probe  will  be  discussed  in  detail 
in  Chapter  IV.  The  cones  on  the  side  walls  are  eight 
inches  from  base  to  tip,  constructed  of  a  carbon 
impregnated  euretnane  material,  and  mounted  on  plywood 
with  contact  cement.  The  first  horizontal  section  of 
the  left  side  wall  (6'  x  4')  was  designed  and 
constructed  with  guides  on  the  top  and  bottom  which 
allowed  forward  and  back  movement  or  entire  removal  of 
the  plywood/cone  side  wall.  The  chamber  was  on  two 
bakelite  covered  work  benches.  The  aluminum  ground 
plane  was  elevated  from  the  workbench  base  with  2x4 
boards  to  provide  a  space  for  probe  telemetry  as  well 
as  coaxial  cable  and  couplings  for  the  antenna 
mounted  vertically  in  the  ground  plane.  Leveling  of 

the  ground  plane  was  accomplished  with  shims  and 
adjustment  of  the  support  legs. 


VARIABLE 


The  following  discussion  includes  the  equip¬ 
ment  involved  in  the  propagation  of  the  EM  wave, 
detection  of  the  IR,  and  the  particular  targets  of 
interest  in  this  report. 

Microwave  Radiation 

Continuous  wave  sources  at  frequencies  of 
937  MHz  and  2.45  GHz  were  used.  The  experimental 
arrangement  is  shown  in  a  block  diagram  in  Figure  48 
and  a  photograph  in  Figure  49.  Table  1  provides 
further  information  on  the  equipment.  As  is  standard 
practice,  tuning  was  accomplished  so  as  to  maximize 
transmitted  power  and  minimize  reflected  power. 

Maximum  source  power  was  80  watts  and  120  watts  for 
the  937  MHz  and  2.45  GHz  sources  respectively.  A 
Hewlett  Packard  spectrum  analyzer  (Model  8555A)  was 
used  to  confirm  that  the  output  frequencies  were  937 
MHz  and  2.45  GHz  for  the  respective  microwave  sources. 
The  reasons  for  selecting  937  MHz  included  remaining 
down  slightly  from  the  limit  of  the  device,  it  had 
been  used  at  Monterey,  and  it  provided  a  significantly 
different  wavelength  than  that  of  the  2.45  GHz  single 
frequency  source.  In  Figure  50  is  shown  a  variable 
aperature,  corner  reflector,  and  a  parabolic  cylinder 
reflector  depicted  with  a  quarter-wave  monopole. 
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Table  1 

Microwave  Equipment 


Nomenclature 
Microwave  Sources 


Isolator 

Dual  Directional 
Coupler 

Attenuators 

Dual  Stub  Tuners 

Power  Meters 


Identification 


Microwave  Power  Generator 
Model  MPG  4M 

KIVA  Instrument  Corporation 
Rockville,  MD  (2.45  GHz) 

Power  Signal  Source 
Sierra  Electronics  Division 
Model  476A-1000 
Philco  Corporation,  USA 
(.5-1.0  GHz) 

Ferrite  Isolator 
Micro-mega 
Model  RL-146 

Hewlett  Packard 
Model  765D 
(450-945  MHz) 

Variable  Attenuators 
Microlab/ FXR 
Model  No.  AJ-310N 
(.5-6.0  GHz) 

Coaxial  Attenuators 
Merr imac 
Type  AU10A 
(.5-12  GHz) 

Microlab 
S  2-05  N 

(Used  at  937  MHz) 

Microlab/FXR 

Model  N  800  A 

Hewlett  Packard 
Model  432A 
Thermister  Mount 
Model  478A 


Coaxial  Cable 


Plastoid  Corporation 
RG-214U 
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The  driven  elements  were  quarter -wave 
monopoles  mounted  on  a  ground  plane.  The  current 
distribution  on  the  antenna  is  modeled  through  image 
theory  as  that  of  a  center  fed  half  wave  dipole  in 
free  space  with  one  half  of  the  radiation  resistance. 
The  antennas  were  made  from  .32  cm  diameter  brass  rods 
and  threaded  so  as  to  be  easily  interchangeable  and 
to  fit  into  a  connector  mounted  into  the  ground  plane 
from  below.  The  interior  line  in  the  coaxial  cable 
was  connected  to  the  antenna  and  maintained  in  a  rigid 
position  with  a  Teflon  spacer  which  fit  firmly  in 
mounted  opening  in  the  ground  plane.  The  accepted 
value  of  .95  was  applied  to  the  electrical  quarter 
wave  length  to  yield  an  antenna  height  of  7.6  cm  and 
2.9  cm  above  the  ground  plane  for  frequencies  of 
937  MHz  and  2.45  GHz  respectively  (29  )•  With 
vertically  oriented  antennas  on  a  horizontal  ground 
plane,  it  is  clear  that  the  incident  field  at  a 
target  would  be  rearly  vertically  polarized,  S  vertical 
and  H  horizontal. 

To  increase  power  at  the  target  and  thus 
increase  perceptible  heating  effects,  greater  gain 
from  the  antenna  was  required.  Additionally,  a  plane 
wave  was  also  required  if  correlations  were  to  be 
reasonably  made  with  theoretical  studies  whicn 


81 


assumed  an  incident  plane  wave.  Corner  reflectors 
and  parabolic  cylinder  reflectors  were  logically 
choices.  It  was  also  recognized  that  a  chamber 
width  of  only  four  feet  could  result  in  reflection 
problems  and  thus  difficulty  in  achieving  a  "quiet 
zone"  at  the  target.  A  hinged,  variable  aperature 
corner  reflector  was  designed  and  constructed  using 
emperically  determined  optimum  dimensions  from  work 
by  Cottony  and  Wilson  (14).  Antenna  to  apex  spacing 
consideration  were  originally  discussed  in  1940  by 
Krauss  (29).  It  is  well  known  that  corner  reflection 
with  aperature  angles  such  as  180°,  90°,  60°,  45° 
etc.,  that  result  when  360°  is  evenly  divided  by  an 
integer  can  be  modeled  through  image  theory  to  readily 
give  the  far  field  pattern  of  interest.  Also,  apex 
to  antenna  spacing  can  be  varied  to  give  a  maximum 
lobe  on  the  centerline,  a  null  on  centerline  or  a 
transition  situation  between  these  extremes  (27). 
However,  the  problem  remained  essentially  an  emperical 
one  due  to  the  complicating  "small"  anechoic  chamber, 
particularly  in  the  case  of  937  MHz.  A  ripple  tank 
study  of  two  dimensional  modeling  of  wave  inter¬ 
actions,  varying  reflectors  and  apex  to  antennas  spac¬ 
ing  in  a  "small"  chamber  was  conducted  and  a  summary 
appears  in  Appendix  A. 


The  far  field  or  Fraunhoffer  condition 


requires  that  the  distance  of  the  scattering  object 
from  an  exciting  source  (D)  be  related  as 


where 

X  is  the  wavelength  and 

L  is  the  characteristic  length  which  in  the 
case  of  a  corner  reflector  is  simply  twice  the  apex 
to  antenna  distance  (See  Figure  51) .  A  90°  corner 
reflector  was  eventually  used  for  the  937  MHz  case 
(X  =  32  cm)  such  that 

S  =  . 32X 

which  implies 

L  =  2S  =  .  64  X  (47) 

and  therefore 


0  >  2 ( . 64 X ) 2 
X 


.  82x 


D  >  26.2  cm 

The  target  was  actually  in  the  chamber  100.5  cm  from 
the  antenna  with  S  =  10.2  cm  implying  that  D  =  110.7  cm 
which  is  in  the  far  field. 

The  positioning  of  a  parabolic  cylinder 
reflector  for  the  2.45  GHz  case  was  again  impir ically 
arrived  at  after  having  started  with  situation  where 


the  antenna  was  at  the  focal  point  and  the  antenna 
was  roughly  in  the  plane  of  the  aperature.  Final 
antenna  to  reflector  spacing  was  15.3  cm.  Free  field 
plots  for  both  frequencies  across  20  cm  at  the  target 
location  for  3/4  and  1  3/4  inches  off  the  ground  plane 
are  found  in  Chapter  V. 

Infrared  Detection 

This  work  extends  efforts  performed  at  the 
Naval  Postgraduate  School  under  Dr,  Robert  W.  Burton 
whore  qualitative  studies  were  carried  out  on 
non-metallic  surfaces. 

Thermography  Equipment 

The  system  used  in  all  of  our  measurements  was 
an  AGA  680  camera  operating  in  the  2-5.6  pm  spectral 
region.  This  device  has  an  8°  lens  with  a  1.3  mr 
instantaneous  field  of  view.  The  detector  is  made  of 
indium  antimonide  (InSb)  and  cooled  to  77  K  with 
liquid  nitrogen  with  a  limiting  sensitivity  of  ,1°K. 
Since  the  detector  is  a  semi-conductor,  it  is  a 
photon  detector  rather  than  an  energy  detector.  This 
is  significant  since  the  number  of  photons  a  blackbody 
emits  is  proportional  to  T^  rather  than  T4  as  is  the 


case  for  the  total  energy  radiated.  The  AGA  680  is  a 
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line  scanning  system  that  forms  a  128  x  128  array  of 
its  8°  square  window.  One  complete  scan  or  "frame" 
takes  1/16  second  and  is  comprised  of  64  rows  having 
128  elements  in  each  row.  It  detects  the  number  of 
photons  emitted  from  each  element  per  unit  time  and 
then  assigns  it  a  level  from  1  to  10  based  on  the 
number  of  photons  received.  These  levels  are  graph¬ 
ically  displayed  in  the  form  of  ten  discrete  colors 
on  a  CRT  display.  The  system's  sensitivity  is  dis¬ 
cretely  variable  in  the  form  of  fixed  isotherm  unit 
(ISU)  windows  which,  in  other  words,  determine  the 
temperature  "band  width"  that  the  camera  will  be 
capable  of  detecting.  It  also  contains  a  threshold 
adjustment  which  determines  where  on  the  temperature 
scale  the  temperature  window  will  be  located.  In 
Figure  52  the  AGA  680  Thermovision  System  is  seen  with 
the  black  and  white  monitor  to  the  left  of  an 
Hewlett  Packard  9845B  minicomputer,  and  the  color 
monitor  to  the  left  of  the  black  and  white  monitor. 

The  camera  is  in  the  foreground  and  was  adequate  for 
this  phase  of  the  overall  project  but  would  not 
provide  the  flexibility  of  placement  near  a  target 
object  due  to  its  size  and  field  perturbing  effects  in 
chamber.  Initial  work  was  required  on  the  system  in 
the  spring  and  summer  of  1980  as  a  faulty  video 
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amplification  board  and  "chopper"  adjustments  were 
located  and  corrected.  However,  even  with  a  properly 
functioning  680  system  the  results  were  still  in 
gray  levels  on  the  black  and  white  display  and  ten 
color  levels  on  the  color  display.  As  will  become 
clear,  a  photographic  display  was  not  adequate  to 
explore  the  surface  current  problem. 

Digitizing  the  Thermovision  Output 

For  the  required  quantitative  analysis,  such 
as  the  point-by-point  correction  for  emissivity 
variations,  and  determination  of  surface  current 
densities,  computer  support  was  necessary.  Unfortu¬ 
nately,  the  data  from  the  Thermovision  camera  was  too 
fast  for  direct  input  to  a  reasonably  compact  computer. 
The  required  interface  between  the  camera  and  the 
computer  was  designed  specifically  for  the  present 
needs  of  the  project.  The  interface  stores  one  frame 
of  data  at  camera  speeds  and  then  off-loads  the  data 
at  computer  speeds  (interface  built  by  Major  Warmuth) . 

The  AGA  Thermovision  680  camera  scans  from 
left  to  right  at  a  rate  of  0.625  ysec  per  line.  It 
uses  approximately  64  lines  (from  top  to  bottom)  to 
produce  one  frame,  and  interlaces  seven  frames  per 
picture.  The  first  frame  illuminates  dots  in  lines 
one,  three,  five,  etc.  The  second  frame  illuminates 


dots  in  lines  two,  four,  six,  etc.  Tne  TV  uses  the 
standard  electrical  power  frequency  to  control 
vertical  movement.  Therefore,  it  produces  sixty 
frames  and  thirty  pictures  per  second.  The  hermo- 
vision  produces  about  16  frames  and  2.3  pictures  per 
second . 

The  interface  reads  128  samples  per  scan 
line,  digitizes  the  temperature  level,  and  stores  the 
results  for  one  frame.  This  yields  approximately  8200 
pieces  of  data  in  62  msec.  The  data  is  then  trans¬ 
ferred  into  a  Hewlett  Packard  9845B  desk-top  computer. 
Once  the  data  is  in  the  computer,  it  may  be 
manipulated  through  normal  matrix  techniques. 

The  computer  signals  the  interface  when  a 
frame  of  data  is  to  be  taken;  the  interface  then  waits 
for  a  vertical  snyc  pulse  from  the  camera.  The 
vertical  sync  pulse  aligns  all  circuitry  in  the  Thermo¬ 
vision  system  at  the  upper  left-hand  corner  of  the 
display.  When  the  interface  receives  the  vertical 
sync,  it  acknowledges  the  computer's  request  and 
begins  digitizing  and  storing  the  video  data.  A  very 
fast  16K  by  8  bit  memory  board  was  locally  available; 
therefore,  design  of  a  new  storage  device  was  not 
required.  Digitizing  initially  to  ten  levels  was 


appropriate  in  the  original  system  since  the  Thermo¬ 
vision  color  display  depicts  only  ten  levels. 

The  interface  uses  an  onboard  clock  running 
at  277  KHz.  This  clock  provides  173  samples  per 
line,  but  about  25%  of  the  line  is  a  right-edge  bar 
and  horizontal  blanking  which  provides  no  useful 
information.  Thus,  only  the  first  128  words  are 
stored.  The  counter  which  controls  word  storage  is 
reset  at  128  and  enabled  by  a  horizontal  sync  pulse 
generated  by  the  Thermovision  camera.  The  horizontal 
sync  pulse  aligns  the  Thermovision  system  at  the  left 
hand  edge  of  the  display.  The  output  of  the  counter 
is  used  to  determine  the  least  significant  eight 
bits  of  the  address  to  the  memory  where  the  current 
digitized  value  of  the  video  signal  is  to  be  stored. 
The  most  significant  six  bits  of  the  address  are 
determined  by  counting  the  number  of  horizontal  sync 
pulses  since  the  last  vertical  sync  pulse.  The  inter 
face  continues  digitizing  and  storing  data  until  it 
sees  a  second  vertical  sync  pulse. 

The  second  vertical  sync  pulse  shuts  down  the 
data  collection  portion  of  the  interface  and  forces 
the  system  to  notify  the  computer  that  a  frame  of 
data  has  been  captured.  When  the  computer  receives 


this  notification,  it  takes  control  of  the  interface 
and  the  associated  memory.  The  computer  transfers  the 
first  address  to  the  memory,  requests  the  data  stored 
there,  and  places  it  in  its  internal  memory.  After 
restoring  the  data,  it  goes  to  the  second  address,  and 
so  on.  When  the  entire  frame  of  data  has  been  moved 
to  the  internal  memory,  the  operator  is  notified  so 
data  manipulation  may  begin,  rot  data  manipulation 
considerations  and  the  resolut ir  constraint  of  taking 
only  one  frame  at  a  time,  the  128  words  per  row  were 
reduced  to  64  by  taking  only  every  other  one.  There¬ 
fore,  the  working  matrix  was  64  x  64. 

In  Figures  53  and  54  are  shown  two  early  efforts 
involving  the  microwave  illumination  of  square  plates 
and  the  digitization  of  the  results.  The  square  plate 
was  a  boron-epoxy  substrate  with  a  carbon-based  paint 
for  a  coating.  The  frequency  was  changed  from  937  MHz 
to  2.45  GHz  such  that  the  electrical  size  of  plate 
changed  from  .50  to  1.3  wavelength.  The  accompanying 
patterns  also  changed  as  is  clear  from  the  figures. 
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Thermovision  Calibration 

With  Thermovision  data  in  digital  form  the 
next  task  is  to  transform  this  data  in  isotherm  units 
to  that  of  temperature  which  will  later  be  shown  to 
be  related  to  current.  Calibration  curves  for  the 
AGA  Thermovision  680  may  be  modeled  by 


i  =  -£■ 


where 


<  48  ) 


I  =  Isotherm  unit  corresponding  to  T 
T  =  Absolute  temperature  (°K) 
p, q  =  Aperature  constaints 
A  typical  calibration  taken  from  680  manual  is  given 
in  Figure  55.  By  measuring  two  values  on  the  curve 
(i.e.,  Tjl /  1^  and  T2,  I2)  and  substituting 
individually  in  equation  (48),  the  result  is  a  system 
of  two  equations  in  two  unknowns.  Thus  taking  two 
measurements  using  a  suitable  temperature  source  can 
yield  the  calibration  unknowns  p  and  q. 

From  equation  <4G}  and  the  two  measured  values, 
the  following  expressions  are  obtained. 


_9L 

P  =  Ij_  (e  Ti  -1) 
q  =  T2  In  +  1) 


(  49  ) 


(  50) 


AGA  THERMOVISION^  SYSTEM  680 
Typical'callbration  chart  -30°  to  + 190° C 


temperature  in  *c 

Figure  55:  Published  calibration  curve  for  the 

Thermovision  680  (1ST*  vs  Temp). 


Substituting  equation  (49)  for  p  into  equation  (50) 
for  q,  the  result  is  the  following  transcendental 
equation . 


q 


T2  in 


(e 


S. 

T1 


-1) 


+1 
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Typically,  q  is  much  larger  than  T^  so  that 


_SL 

(e  T1  -1) 


and  therefore 


q 


T1T2 

VT2 
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An  example  of  the  relative  magnitude  of  q  and  is 
found  by  taking  two  points  from  the  graph  in  Figure 
55,  Tx  =  90°  C  =  363°  K,  I1  =  150  ISU,  T2  =  190°  C  = 
463°  K,  and  I2  =  950°  and  solving  for  p  and  q.  The 
resulting  expression  defining  the  curve  is 


I 


771,348 

3102 


-1 


where 


p  =  771,348  and 
q  =  3102 

A  more  general  solution  to  the  transcendental 
equation  is  found  in  computer  program  named  'CAMCAL' 


in  Appendix  D  using  the  Newton's  approximation.  This 


itterative  technique  allows  the  user  to  select  the 
convergence  criterion  and  uses  the  following 


function. 

f(q) 


T2  In 


+1 


-q 


its  derivative 


f'(q)  = 


and  the  standard  recursion  relation. 


(54; 


(55; 


dn+l 


-<n-£!V 


£,<V 
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The  values  of  p  and  q  used  in  all  calculations  in 
this  paper  were  determined  to  be 


p  =  772600.6  and 
q  =  3102.6 

which  yield  the  calibration  curve  in  Figure  56.  Note 
that  for  ambient  temperature  and  for  small  changes  in 
temperature*  the  relationship  between  isotherm  units 
and  degrees  of  temperature  is  nearly  linear. 


Targets 


Several  considerations  were  involved  in  the 
selection  of  flat  plate  electrical  and  physical  sizes. 
The  minimum  focal  distance  for  the  AGA  Thermovision  680 
is  1.7  m  wnich  was  used  in  data  collection.  Thus  the 
object  size  was  reasonably  restricted  to  greater  than 
10  cm  along  the  largest  physical  dimension  to  obtain  a 
sufficient  number  IR  elements  to  distinguish  current 
patterns.  For  the  electrical  dimension,  chamber  size 
considerations  reasonably  restricted  a  "quiet  zone"  to 
approximately  20  cm  for  the  maximum  horizontal  direc¬ 
tion.  One  must  recall  that  objects  in  contact  with  a 
ground  plane  will  have  an  electrical  image  such  that 
an  electrical  aguaie  plate  will  physically  appear  as 
a  rectangular  plate  having  a  horizontal  dimension  twice 
that  of  the  vertical  dimension.  Additionally,  the 
investigation  of  current  patterns  should  include 
resonant  and  nonresonant  sizes.  The  following 
physical  and  electrical  sizes  were  chosen. 

Table  2 

Flat  Plate  Size 

Largest  physical  2.45  GHz  937  MHz 

dimension  (cm) _ (Wavelengths ) _ (Wavelengths ) 

16  1.333  .50 


12.24 


1.0 


.3825 
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Three  material  compositions  were  used.  The 
first  was  1/16  inch  plexiglass  coated  on  one  side  with 
a  carbon  based  paint  known  commercially  as  Television 
Tube  Kote  (aquadac) .  The  thickness  cf  several 
sprayed  on  coats  was  59  ym  and  the  conductivity  of 
the  aquadac  was  measured  by  Major  Vic  Martin  to  be 
315  mhos/m.  To  achieve  a  good  electrical  contact  with 
the  ground  plane,  the  plates  were  first  mounted  on 
square  1/32  inch  thick  aluminum  sheet  with  silver 
conducting  point  as  seen  in  Figure  57.  When  placed  in 
the  chamber  the  mounting  aluminum  sheet  was  taped  to 
the  ground  plane  with  copper  tape.  The  second 
composition  was  again  developed  by  Martin.  A  2  mm 
thick  coating  of  carbon/paraffin  mixture  with  a 
measured  conductivity  of  1  mho/m  was  placed  on  a  1/2 
inch  thick  sheet  of  styrofoam  which  in  turn  was 
attached  with  double-stick  material  to  a  1/16  inch 
t'.xck  aluminum  plate.  A  sampling  of  this  composition 
is  seen  in  Figure  58.  The  third  type  was  simply  1/4 
inch  thick  aluminum  plates  that  were  sanded  to  remove 
a  factory  coating  and  used  for  the  probe  measurement 
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(b)  Squace  flat  plate. 

Figure  57:  Targets  -  Aquadac  on  plexiglass 


CHAPTER  IV 


MAGNETIC  FIELD  PROBE 

The  classic  technique  for  measuring  surface 
currents  is  with  a  loop  antenna  (50) .  In  most  cases, 
the  objects  are  metallic  with  a  structure  which  is 
convenient  for  shielding  the  probe  telemetry  from  the 
incident  and  scattered  EM  field.  it  was  not  clear 
during  the  work  on  this  project  that  the  taraet  would 
be  restricted  to  metal.  A  magnetic  field  probe  was 
therefore  developed  that  could  scan  along  any 
material  surface  or  be  used  as  a  free  field  probe. 

The  basis  of  operation  of  a  magnetic  field  or 
"current"  probe  is  Faraday's  Law.  Stated  roughly, 
the  induced  voltage  in  the  loop  will  be  equal  to  the 
time  rate  of  change  of  the  magnetic  field  through  the 
loop  assuming  permeability  and  area  of  the  loop  are 
constant.  Most  exactly,  the  line  integral  of  the 
electric  field  around  the  perimeter  of  a  loop  is 
equal  to  the  integral  of  the  time  rate-of-change  of 
the  magnetic  induction  over  the  area  of  the  loop. 
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where  E  is  the  electric  field  intensity  in  volts  per 
meter#  and  B  is  the  magnetic  flux  density  in  teslas. 
The  left  side  of  equation  is  simply  the  induced 
voltage  (V)  in  the  loop.  The  right  side  can  be 
simplified  by  assuming  that  B  is  constant  over  the 
area  S.  Assuming  a  sinusoid  field,  the  greatest 
dimension  for  any  probe  in  this  report  is  .64  cm 
compared  with  the  smallest  wavelength  of  12.24  cm 

t  <  — V 

implies  that  the  variation  of  H  over  the  area  will  be 
no  larger  than  five  percent.  Therefore,  the  component 
of  6  parallel  to  d3  can  be  considered  spatially 
independent  of  d§  over  the  loop  and  the  area  is 
independent  of  time  such  that  equation  (57)  reduces 
to 


V  =  '  S  Hi  (58) 

where  B  now  is  the  component  perpendicular  to  the 
plane  of  the  loop.  Applying  the  constituent  relation 
B  =  pH  to  the  case  of  the  probe  in  free  space  (u  =  uQ  = 
permeability  of  free  space) ,  we  have 


„  G  dH 

V  =  -  y  S  -rr 
o  dt 


(59) 


or 


V  =  -  jwy  SH 


(60) 


for  an  e^wt  time  variation  of  the  magnetic  field. 
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Vrms  actually  the  measured  quantity.  Now,  a 
magnetic  field  prohe  can  be  used  to  measure  surface 
current  density  since  the  magnetic  field  over  a 
conductive  sheet  is  related  to  the  surface  current 
density  by  Ampere’s  Law: 


^H*dt  =  /jt*ds  , 


(61) 


where 

H  =  magnetic  field  intensity  in  amperes 
per  meter 

dfl,  =  incremental  length  on  a  closed  contour 
in  meters 

Jt  =  conduction  current  density  plus  the 

displacement  current  density  in  amperes 
per  meter 

ds  =  incremental  area  on  the  contour  in 
square  meters. 

For  the  case  of  metallic  sheets,  the  conduction 
current  density  dominates  the  displacement  current 
density  and  the  magnetic  flux  does  not  appreciaoly 
penetrate  the  sheet.  Therefore,  Ampere's  Law 
reduces  to 

H  =  Js  (62) 

where  Js  is  the  surface  current  density  in  amperes  per 


meter . 


For  metal  plates,  the  assumptions  associated  with  a 
small  probe  and  a  good  conductor  allow  one  to  see 
quite  simply  that  output  voltage  from  a  magnetic 
field  sensing  probe  positioned  at  the  surface  of  a 
conductor  is  linear  with  surface  curren  density. 

If  there  is  appreciable  penetration  of  the 
magnetic  field  through  the  object  surface  the 
simplifying  assumptions  should  be  checked  for 
validity.  In  this  report,  the  probe  will  only  be 
used  on  a  high  conductive  aluminum  surface  in  a  high 
frequency  (gigahertz)  range  which  implies  a  very 
small  skin  depth  (6)  which  is  the  depth  in  the 
conductive  sheet  where  the  current  density  is  1/e  = 
.368  of  its  value  at  the  surface.  The  skin  depth  is 
given  by 

s  *  O'*  1631 

An  important  property  of  a  field  probe  is 
that  it  is  an  integrating  device.  That  is,  the 
contributing  current  density  source  for  the  magnetic 
field  and,  hence  the  induced  voltage,  is  not  a 
"point"  source,  but  rather  an  area  in  the  vicinity  of 
the  probe.  A  generally  accepted  figure  for  probe 
integration  area  is  1.5  times  the  loop  area.  This 
integrating  property  is  particularly  evident  at 
plate  edges  where  it  should  be  understood  that  the 
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current  sources  are  not  present  throughout  the 
integrating  area  of  the  probe  as  would  be  the  case 
in  the  center  of  the  plate.  Experimental  values  of 
current  will  be  diminished  somewhat  from  the 
anticipated  higher  edge  current  values  due 
predominantly  to  this  effect. 

Experimental  Considerations 

Receiving  and  recording  the  induced  voltage 
in  a  high  frequency  (EM)  field  presents  many 
experimental  problems.  The  only  excitation  source 
that  is  desired  is  the  time  varying  magnetic  field 
at  the  probe.  The  potential  coupling  of  the  electric 
field  to  the  probe  and  potential  coupling  of  both 
£  and  H  to  the  telemetry  must  be  overcome.  The 
recording  device  must  be  capable  of  either  receiving 
at  the  driving  frequency  or  the  signal  must  be 
converted  to  D.C.  The  strength  of  the  incident 
field  and  the  sophistication  of  equipment  dictated 
the  minimum  reasonable  probe  size  and  therefore  the 
minimum  integrating  area. 

Three  probes  were  developed  with  varying 
degrees  of  success  and  are  shown  in  the  following 
three  figures.  The  AC  semirigid  coaxial  approach 


(Figure  59)  utilized  the  cable  from  Uniform  Tubes 


Figure  59:  Photograph  of  a  Magnetic  Field 

Probe  -  Semirigid  Coaxial  Cable 
Construction . 
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Figure  60:  Drawing  of  a  Magnetic  Field 

Probe  -  Semirigid  Coaxial 
Cable  Construction. 


Figure  63: 


Photograph  of  the  Horizontal 
Magnetic  Field  Probe. 
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Figure  64: 


Drawing  of  the  Horizontal 
Magnetic  Field  Probe. 
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Company,  which  was  connected  to  an  RG  58  C/U  cable 
and  finally  into  an  HP  8555  A  Spectrum  Analyzer 
RF  Section  with  an  HP  855  A  Spectrum  Analyzer  - 
IF  Section  and  141T  Display  Section.  The  D.C. 
approach  ^Figures  61  and  63) ,  used  a  Schottky-barrier 
diode  as  a  portion  of  the  loop  with  the  signal 
carried  by  a  braided  shield,  coaxial  cable  with  a 
measured  capacitance  of  87  pF/ft.  Horizontal  and 
vertical  orientations  of  the  probes  corresponded  to 
horizontal  and  vertical  current  components.  The  AC, 
semirigid  coaxial  approach  involved  more  difficulty 
with  mechanical  scanning,  required  an  operator  to 
read  and  record  from  a  CRT,  and  was  only  used  as 
rough,  independent  verification  of  the  D.C.  results. 
Fabrication  of  the  loop/diode  system  evolved  to  a 
semi-automatic  scanning  procedure.  The  orientation 
of  the  probe  and  coaxial  cable  in  the  chamber  is 
seen  in  Figure  65.  Recall  that  the  incident  EM 
field  is  vertically  polarized  such  that  a  horizontal 
cable  has  little  effect  on  the  wave.  For  structural 
and  grounding  considerations  the  outside  shielding 
on  left  and  right  side  of  the  loop,  as  viewed  down  the 
chamber,  were  soldered  together.  The  loop  was 
designed  with  the  diode  positioned  at  the  target 
surface  where  the  boundary  conditions  imply  a  minimum 
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horizontal  component  of  the  electric  field.  The 
loop  was  formed  from  the  diode  wire,  covered  with 
insulating  shriek-tape  up  to  the  diode,  and 
soldered  to  the  center  wire  of  the  coaxial  cable. 
Except  for  the  diode  casing,  the  loop  wire  and  center 
coaxial  cable  lead  were  coated  with  epoxy  cement, 
let  dry,  then  painted  with  silver  conducting  paint 
which  also  made  a  good  electrical  contact  with  the 
outer  coaxial  cable  shielding.  The  signal  was 
rectified  by  the  capacitance  of  the  coaxial  cable 
before  being  amplified. 

The  recording  equipment  for  the  magnetic 
field  probe  signal  is  seen  in  Figure  66.  On  top  of 
the  cart  is  a  Hewlett  Packard  (Model  7004B)  X~Y 
recorder  with  D.c.  preamplifiers  (Model  17171A)  for 
both  X  and  Y  inputs.  Prior  to  the  signal  reaching 
the  recorde-,  it  was  sent  through  an  amplifying 
circuit  located  in  the  metal  box  under  the  recorder 
which  provided  shielding.  A  common  ground  for  the 
coaxial  cable  snield,  the  Proto-Board  containing 
the  circuit,  the  shielding  metal  box,  and  the 
recorder  was  provided  with  10  AWG  copper  wire 
attached  to  a  water  pipe.  The  reasons  for  building 
the  circuit  shown  in  Figure  67,  were  for  the  high 
impedance  (1.5  TO),  provided  by  the  operational 


amplifier  (CA  3140  -  FET  device),  additional  filtering 
for  60  cycle  noise,  and  the  gain  provided  prior  to 
reaching  the  Y  input  of  the  X-Y  recorder  by  way  of 
another  shielded  cable. 

The  X-input  to  the  X-Y  recorder  provided 
horizontal  position  information  about  the  probe  in 
the  chamber.  In  Figure  68  is  seen  a  close-up  of  the 


Figure  68:  Experimental  arrangement 

for  horizonLal  and  vertical 
probe  positioning. 
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pulley  assembly  showing  the  coaxial  cable  existing 
through  a  slit  cut  in  the  side  wall  of  the  chamber. 
After  passing  over  the  upper  pulley,  which  is 
matched  on  the  opposite  side  of  the  chamber,  the 
cable  passes  over  a  rubber  covered  pulley  attached 
to  the  shaft  of  a  IK  resistance,  10  turn  potentio¬ 
meter  (Helipot) .  With  a  5V  power  supply  connected 
across  the  potentiometer,  and  a  vernier  feature  on 
the  X-input  pre-amplifier  on  the  recorder,  one-to- 
one  distance  information  was  provided  between  probe 
position  and  X-coordinate. 

The  mechanical  scanning  system  was  simple 
and  effective.  The  horizontal  probe  movement  was 
accomplished  using  the  1/40  HP  variable  speed, 
reversible  laboratory  stirrer  (G.K.  Heller  Corpora¬ 
tion,  Model  GT-18  connected  to  a  GT-21  Motor 
Controller)  shown  connected  to  a  lead  plate  in  Figure 
66.  A  casting  line  was  connected  to  a  strap  mounted 
on  the  coaxial  cable  and  to  a  threaded  shaft  on  the 
motor  which,  when  counter  balanced  from  the  opposite 
side  of  the  chamber,  provided  for  positive,  variable 
speed  horizontal  displacement  of  the  probe.  The 
vertical  position  was  simply  adjusted  by  placing 
matched,  horseshoe  shaped,  wooden  spacers  under  the 
vertical  stops  of  the  pulley  as  seen  in  Figure  68. 
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Calibration 

Several  measurements  were  taken  to  establish 
the  linearity  of  the  recorded  induced  voltage  to  the 
magnetic  field  strength  at  the  probe.  The  first 
step  was  to  vary  a  D.C.  input  voltage  to  the 
magnetic  field  probe  pre-amplifier  circuit  and 
record  the  output.  The  results  are  shown  in  Figure 
69  (a  gain  of  5  was  eventually  sufficient).  The 
next  step  was  to  measure  the  power  at  the  probe  and 
compare  it  to  the  recorded  voltage.  Two  calibrated 
power  density  meters  (Figure  70)  were  available  and 
tested  against  each  other  (a  General  Microwave 
Corporation  Radiation  Hazard  Meter,  Model  481A 
with  probe  Model  82  and  a  NARDA  Microline  Electro¬ 
magnetic  Monitor,  Model  8100  with  probe  Models  8122A 
and  8121A) .  The  situation  of  needing  a  power  meter 
in  the  chamber  while  taking  probe  measurements  was 
circumvented  by  establishing  linearity  of  the 
microwave  source  output  meter  readings  with  that  of 
a  power  meter  located  at  the  target  location  with 
both  target  and  probe  removed.  The  results  are 
shown  in  Figures  71  through  74 .  Measurements  were 
now  taken  of  the  X-Y  recorder  response  to  varying 
microwave  source  power  inputs.  The  data  best  fit 
with  a  second  degree  polynomial  (Figures  75  to  78) 


MEASURED 


Figure  70: 


Microwave  power  meters. 


as  one  would  expect  since  power  is  related  to  the 
magnetic  field  as  follows: 


(E  x  H*) 


H2 


since 


E 

H 


in  free  space 


S  =  time  averaged  Poynting  vector 


where 


$  =  electric  field  intensity 

H  =  magnetic  field  intensity 
Vo  =  peameability  of  free  space 
c0  =  permittivity  of  free  space 

The  relationship  of  the  square  root  of  power  to 
induced  voltage  at  the  recorder  is  seen  in  Figures 
79  through  82  to  be  linear  in  the  region  of 
experimental  interest  (nonlinear  diode  characterstics 
are  evident  very  near  zero) . 

Very  similar  results  were  found  using  the 
following,  more  sophisticated  equipment  from  the 
Precision  Measurement  Equipment  Laboratory: 

Power  Meters  -  HP  432A  (8478B  Thermistor 

Mount) 

Signal  Generator  -  HP  8640B  (937  MHz) 

-  HP  8620C  Sweeper  Oscillator 
with  86222B-H69  RF  plug-in 
The  signal  source  was  connected  at  the  diode  and 
the  output  from  a  microvoltmeter  was  observed  at 
the  X-Y  recorder.  Results  appear  in  Figures  83 
through  90. 

Therefore,  the  voltage  seen  at  the  X-Y 
recorder  is  linear  with  the  magnetic  field  at  the 
probe  and  thus,  linear  with  surface  current  density 
on  the  metal  plates. 
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Figure  71:  Correlation  of  microwave 

source  power  indication 


to  free  field  power  measured 
at  target  location  (937  MHz) . 


Figure  72: 


Plot  of  source  power  vs  free 
field  power  (937  MHz). 
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probe  response  to  source 
power  (937  MHz) . 
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Figure  77:  Correlation  of  magnetic  field 

probe  response  to  source 
power  (2.45  GHz) . 
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Figure  78:  Plot  of  magnetic  field 

probe  response  vs  source 
power  (2.45  GHz) . 
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power  (2.45  GHz) . 


Figure  81 


Figure  82:  Plot  of  magnetic  field 

probe  response  vs  square 
root  of  source  power 
(2.45  GHz) . 
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Correlation  of  magnetic  field 
probe  response  to  input  power 
(937  MHz) . 


Figure  84:  Plot  of  magnetic  field 

probe  response  vs  input 
power  (937  MHz) . 
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Figure  85:  Correlation  of  magnetic  field 

probe  response  to  input  power 
(2,45  GHz)  . 
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Figure  87:  Correlation  of  magnetic  field 

probe  response  to  square  root 
of  input  power  (937  MHz) . 


Figure  88: 


Plot  of  magnetic  field 
probe  response  vs  square 
root  of  input  power  (937 
MHz)  . 
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ment  Equipment  Laboratory  (USAF  Academy) . 


89;  Coorelation  of  magnetic  field 
probe  response  to  square  root 
of  input  power  (2.45  GHz). 
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Figure  go; 

Plot  of  magnetic  field 

• 

probe  response  vs  square 
root  of  input  power  (2.45  GHz). 

* 

• 

-  25  H 
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CHAPTER  V 

PROBE  MEASUREMENT  OF  SURFACE  CURRENTS 

To  obtain  total  current  at  any  point  on  an 
object,  two  orthogonal  components  are  required.  For 
a  plane  polarized  wave  at  normal  incidence  on  simple 
symmetric  flat  plates,  it  is  well  known  that  certain 
"cuts"  or  scans  across  the  plate  will  have  a 
dominant  current  component  in  the  direction  of  the 
incident  electric  field.  The  terms  "vertical"  and 
"horizontal"  will  continue  to  refer  to  the 
experimental  arrangement  whereby  incident  l!  is 
vertical  and  $  is  horizontal.  The  semirigid  coaxial 
cable  probe  measurements  roughly  corresponded  to  the 
loop/diode  measurements  and  only  the  loop/diode  type 
probs  are  considered  here. 

Procedure 

The  probe  scans  were  done  at  three  horizontal 
positions  on  each  of  the  plates. 

The  height  above  the  ground  plane  were 
consistently  taken  as  follows: 


2 


height  (inches) 

Scan 

height 

(inches) 

Upper 

Middle 

Lower 

2.41  (6.12  cm) 

2 

1.25 

.5 

3.15  (8  cm) 

2.5 

1.5 

.5 

These  heights  were  selected  so  as  to  sample  the 
surface  currents  in  unique  locations  with  respect 
to  probe  integration  areas  and  to  remain  clear  of 
the  needed  copper  tape  used  for  support  and  electric 
contact. 

Precautions  were  taken  to  assure  that  the 
probe  remained  in  a  given  orientation  to  the  plate 

with  an  insulatinq  jacket  in  contact  with  the  plate. 
As  seen  in  Figure  91,  a  styrofoam  support  was  made 
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in  which  the  probe  remained  fixed.  The  plates  were 
positioned  in  an  identical  centerline  location  at 
100.5  cm  from  the  antenna  for  probe  and  IR  measure¬ 
ments.  This  location  was  slightly  ahead  of  the 
coaxial  cable  so  that  a  slight  positive  pressure  was 
present  as  the  styrofoam/probe  structure  was  guided 
along  the  plate  surface. 

Results 

Both  horizontal  and  vertical  probe  measure¬ 
ments  were  made.  The  vertical  orientation  results 
are  presented  in  this  chapter  using  scan  levels 
mentioned  earlier  with  the  exception  of  the  free 
field  scans  which  were  done  at  .75  and  1.75  inches 
from  the  ground  plane. 

The  results  for  the  horizontal  orientation 
are  presented  in  Appendix  E .  Free  field 
samplings  at  .75  and  1.75  were  unable  to  detect  a 
signal  as  would  be  expected  with  the  incident  H  field 
also  in  the  plane  of  the  loop.  The  induced  voltages 
were  down  from  the  dominant,  vertical  orientation  but 
not  quite  as  symmetric,  revealing  possible  unexpected 
coupling. 


ical  Magnetic  Field  Probe  Measurements  for  the 
wavelength  circular  disk. 


937  MHz 
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CHAPTER  VI 

IR  MEASUREMENT  OF  SURFACE  CURRENTS 

The  capability  of  detecting  a  current  induced 
heating  pattern  and  taking  the  digitized  information 
in  isotherms,  converting  to  temperature,  then  finally 
to  currents  will  be  shown.  The  target  compositions 
were  discussed  in  Chapter  III.  It  must  be  emphasized 
that  current  densities  obtained  through  the  IR 
technique  presented  here  are  the  current  densities 
in  the  given  coatings,  either  aquadac  on  plexiglass, 
or  carbon/paraffin  over  metal. 

Technique 

Having  discussed  the  digitization  of  thermo¬ 
vision  data,  and  the  application  of  the  calibration 
curve  to  transform  from  isotherm  to  temperature,  the 
relationship  between  current  density  and  temperature 
must  be  made.  This  is  accomplished  by  fitting 
emperically  obtained  data  to  a  functional  relation,  a 
second  order  polynomial. 

The  experimental  arrangement  for  relating 


current  density  is  given  in  the  following  figure 


layer  of  adquadac  on  a  1/16  inch  thick  plexiglass 
plate.  Two  parallel  copper  tape  strips  were  placed 
on  the  surface  and  a  thin  connecting  strip  of 
silver  paint  was  applied  on  the  copper  and  aquadac 
surfaces  overlapping  the  junction  to  assure  that 
good  electrically  contact  was  made.  The  copper  strip 
electrodes  were  connected  to  a  Kepco  Power  Supply 
(Range  0-36V,  0-1. 5A)  and  a  digital  ammeter  (Sabtronics 
Model  2000).  A  plate  of  identical  composite  was  then 
placed  near  the  test  plate  (right  side  in  Figure  103) 
to  provide  an  ambient  reference  source  of  identical 
emissivity.  Voltage  levels  were  incrementally 
increased  across  the  electrodes  as  the  differential 
temperature  levels  between  test  and  reference  plates 
were  determined  using  the  Thermovision  system.  At 
each  new  total  current  level  the  system  was  allowed 
to  reach  steady  state  (about  15  minutes)  before 
readings  were  taken.  It  was  assumed  that  the  current 
densities  were  uniformly  distributed  throughout  the 
thin  coatings. 

The  data  was  then  fit  to  a  quadratic  function 
and  plotted  using  a  program  called  ’J  vs  T'  found  in 
Appendix  D*  A  maximum  likelihood  routine  was  used 
to  find  the  coefficients  of  the  squared  (Vlavg)  and 
linear  (V2avg)  terms.  The  results  for  the  two  target 


compositions  are  found  in  Figures  104  and  105-  Now, 
on  a  point  by  point  basis,  knowing  the  differential 
temperature  between  an  illuminated  elemental  area 
and  the  ambient  condition  the  functional  relation  was 
easily  solved  for  current  density.  Having  the  needed 
relationship  to  obtain  current  density,  the  problem 
of  random  errors  associated  with  the  interface 
remained. 

A  statistical  procedure  was  developed  to 

reduce  the  influence  of  random  errors  on  the  final 

matrix  from  which  calculations  were  made.  The  IR 

data  was  accomplished  with  known  sources  of  IR, 

such  as  overhead  lights,  extinguished,  and  a  steady 

state  heating  pattern  established.  Approximately  10 

2 

minutes  of  illumination  at  approximately  10  mW/cm 

2 

at  937  MHz  and  approximately  15  mW/cm  at  2.45  GHz 
were  used.  Seven  frames  of  data  were  collected  for 
the  statistical  analysis  in  all  cases.  Sensitivity 
settings  of  one  and  two  were  used  such  that  nearly 
full  scale  variation  of  £  isotherm  units  (0  to  9) 
existed  on  the  surface  but  without  being  outside  of 
the  selected  sensitivity  window.  Any  spurious 
readings  outside  of  the  0  to  9  range  were  not  included 
in  the  averaging  procedure.  The  remaining  data  was 
considered  on  an  element  by  element  basis.  The  values 


at  an  element  position  were  averaged,  a  standard 
deviation  was  taken  and  points  outside  of  one  standard 
deviation  were  discarded.  Thus,  a  new  mean  was 
computed  from  the  remaining  values  and  the  element 
was  now  assigned  this  new  mean  for  the  "FINAL" 
matrix.  This  procedure  is  the  basis  of  the  program 
called  'IRSTAT'  found  in  Appendix  d. 

Results 

From  the  statistical  procedure  the  remaining 
matrix  or  "FINAL"  matrix  could  be  transformed  to  a 
temperature  matrix  or  a  current  matrix  using  the 
previously  discussed  development  which,  along  with 
some  graphics  routines,  comprise  the  basis  of  a 
program  called  "MEGAIR"  found  in  Appendix  D. 

All  electric  sizes  are  presented  for  the 
composition  of  carbon/paraffin  over  metal.  Only  the 
.50  and  1.0  wavelength  disks  and  square  plates  are 
examined  for  the  aquadac  on  plexiglass  case.  Tne 
digital  presentation  in  numbers  for  the  temperature 
and  current  matrices  and  complete  profile  matrices 
for  the  illuminated,  temperature,  and  current  matrices 
are  given  for  the  1.0  wavelength  square  plate  for 
carbon/paraffin  over  metal. 


The  depicted  profiles  are  IR  obtained  matrix 
rows  that  correspond  to  the  horizontal  scans  of  the 
magnetic  field  probe  given  in  Chapter  V.  The  key  is 
as  follows: 

Upper  -  dotted  line  (••?••) 

Middle  -  broken  line  ( — •  •  —  •• — ) 

Lower  -  solid  line  ( - ) 
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Figure  106:  Thermogram  for  an  illuminated 
.38  wavelength  circ  ar  disk  - 
carbon/paraffin  over  metal. 
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Figure  107: 


Digitized  thermovision  data 
for  .38  wavelength  circular 
disk  -  carbon/paraffin  over 
metal . 
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Figure  108:  Thermovision  data  profiles  for 
.38  wavelength  circular  disk  - 
carbon/paraffin  over  metal. 
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Temperature  profiles  for  .38 
wavelength  circular  disk  - 
carbon/paraffin  over  metal. 
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Figure  111:  Thermogram  for  an  illuminated 
.38  wavelength  square  plate  - 
carbon/paraffin  over  metal. 
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Figure  112:  Digitized  thermovision  data  for 
.38  wavelength  square  plate  - 
car  bon/par  a f f in  over  metal. 
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Figure  113:  Thermovision  data  profiles  for 
.38  wavelength  square  plate  - 
carbon/paraffin  over  metal. 
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Figure  115:  Infrared  obtained  current  profiles 
for  .38  wavelength  square  plate  - 
carbon/paraffin  on  metal. 
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Figure  116:  Thermogram  for  an  illuminated 
.50  wavelength  circular  disk  - 
carbon/paraffin  over  metal. 
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Figure  H7;  Digitized  thermovision  data  for 
.50  wavelength  circular  disk.  - 
carbon/par  a f f  in  over  metal. 


10 


■4 


Figure  118:  Thermovision  data  profiles  for  ■' 

.  50  wavelength  circular  disk  - 

carbon/paraffin  over  metal.  3 
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Figure  119:  Temperature  profiles  for  . 

wavelength  circular  disk  - 
carbon/paraffin  over  metal 
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Figure  120:  Infrared  obtained  current 

profiles  for  .50  wavelength 
circular  disk  -  carbon/ 
paraffin  over  metal. 
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Figure  121: 


Thermogram  for  an  illuminated 
.50  wavelength  circular  disk  - 
aquadac  on  plexiglass. 
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Figure  1?2: 


Digitized  thermovision  data  for 
.50  wavelength  circular  disk  - 
aquadac  on  plexiglass. 
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Figure  123:  Thermovision  data  profiles  for 
.50  wavelength  circular  disk  - 
aquadac  on  plexiglass. 
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Figure  125  Infrared  obtained  current 
profiles  for  .50  circular 
disk  -  aquadac  on  plexiglass. 


Figure  126: 


Thermogram  for  an  illuminated 
.50  wavelength  square  plate  - 
carbon/paraffin  over  metal. 
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Figure  127:  Digitized  thermovision  data  for 
.50  wavelength  square  plate  - 
carbon/paraffin  over  metal. 
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Figure  128:  Thermovision  data  profiles 
for  .50  wavelength  square 
plate  -  carbon/paraffin  over 
metal . 
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Figure  130:  Infrared  obtained  current 

profiles  for  ,50  wavelength 
square  plate  -  carbon/ 
paraffin  over  metal. 
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Figure  131:  Thermogram  for  an  illuminated 
.50  wavelength  square  plate  - 
aguadac  on  plexiglass. 
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Figure  133:  Thermovision  data  profiles  for 
.50  wavelength  square  plate  - 
aquadac  on  plexiglass. 
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Figure  135:  Infrared  obtained  current 

profiles  for  .50  wavelength 
square  plate  -  aquadac  on 
plexiglass . 
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Figure  136:  Thermogram  for  an  illuminated 
1.0  wavelength  circular  disk  - 
carbon/paraffin  over  metal. 
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Figure  137:  Digitized  thermovision  data  for 
1.0  wavelength  circular  disk  - 
carbon/paraffin  over  metal. 
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Figure  138:  thermovision  data  profiles  for  1.0 
wavelength  circular  disk  -  carbon/ 
paraffin  on  metal. 
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Figure  140:  Infrared  obtained  current  profiles 
for  1.0  wavelength  circular  disk  - 
carbon/paraffin  on  metal. 
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Figure  141: 


Thermogram  for  an  illuminated 
1.0  wavelength  circular  disk 
aquadac  on  plexiglass. 
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Figure  142: 


Digitized  thermovision  data 
for  1.0  wavelength  circular 
disk  -  aquadac  on  plexiglass 
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Figure  143:  Thermovision  data  profiles 
for  1.0  wavelength  circular 
disk  -  aquadac  on  plexiglass. 
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Figure  146:  Thermogram  for  an  illuminated 
1.0  wavelength  square  plate  - 
carbon/paraffin  over  metal. 
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Figure  147:  Digitized  thermovision  data  for 
1.0  wavelength  square  plate  - 
carbon/paraf fin  over  metal. 
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Figure  148:  Temperature  matrix  for  1.0  wavelength 
square  plate  -  carbon/paraffin  over 
metal . 
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Figure  149:  Current  matrix  for  1.0  wavelength 

square  plate  -  carbon/paraffin  over 
metal . 
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Figure  150:  Thermovision  data  matrix  profile  for 
1.0  wavelength  square  plate  -  carbon/ 
paraffin  over  metal. 
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Figure  151:  Temperature  matrix  profile  for 
1.0  wavelength  square  plate  - 
carbon  paraffin  over  metal. 


Horizontal  <X)  Rxla 


CURRENT  DENSITY  MATRIX  (J) 

Figure  152:  Infrared  obtained  current  matrix 
profile  for  1.0  wavelength  sqi  are 
plate  -  carbon  paraffin  over  metal 
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Figure  155:  Infrared  obtained  current  profiles 
for  1.0  wavelength  square  plate  - 
carbon/paraffin  over  metal. 


Figure  156 


Thermogram  for  an  illuminated 
1.0  wavelength  square  plate  - 
aquadac  on  plexiglass. 
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Figure  157 : 


Digitized  thermovision  data 
for  1.0  wavelength  square 
plate  -  aquadac  on  plexiglass 
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Figure  159:  Temperature  profiles  for  1.0 
wavelength  square  plate  - 
aquadac  on  plexiglass. 
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Figure  161: 


Thermogram  for  an  illuminated 
1.3  wavelength  circular  disk  - 
carbon/paraffin  over  metal. 
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Figure  162:  Digitized  thermovision  data  for 
1.3  wavelength  circular  disk  - 
carbon/paraffin  over  metal. 
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Figure  163:  Thermovision  data  profiles  for 
1.3  wavelength  circular  disk  - 
carbon/paraffin  over  metal. 
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Figure  164:  Temperature  profiles  for 

1.3  wavelength  circular  disk  - 
carbon/paraffin  over  metal. 
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Figure  165:  Infrared  obtained  current  profiles 
for  1.3  wavelength  circular  disk  - 
carbon/paraffin  over  metal. 
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e  169:  Temperature  profiles  for 
1.3  wavelength  square 
plate  -  carbon/paraffin 
over  metal. 
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CHAPTER  VII 


CORRELATION 

Selected  situations  for  relating  current 
density  distributions  froiri  theory,  magnetic  field 
probe  measurements,  and  IR  obtained  results  will  be 
presented.  The  computer  solution  for  the  perfectly 
conducting  disk  was  favorably  compared  in  Chapter  II 
to  disks  of  smaller  and  larger  electrical  size  than 
those  observed  in  the  experimental  portion  of  this 
work.  The  theoretical  solutions  to  the  square  plates 
differ  and  the  only  directly  applicable  scan  to  this 
study  would  be  the  1.0  wavelength  square  plate  and 
the  observation  that  the  small  plate  (.15  wavelength) 
solved  by  Wilton  compares  well  in  form  with  the 
smallest  plate  (.38  wavelength)  used  in  this  work. 

The  .38  wavelength  disk  and  square  plate  provide 
situations  where  the  dominant  current  component  is 
vertical  and  thus  probe  measurement  for  the  vertical 
component  can  be  considered  as  that  of  total  current, 
thus  compared  with  the  IR  result  for  total  current. 


The  disk  problem  has  been  solved  theoreti¬ 
cally  in  Chapter  II  in  the  logical  coordinate  system, 
cylindrical.  To  compare  the  theoretical  disk  results 
with  the  magnetic  field  probe  scans,  the  polynomial 
fits  for  the  maximum  Jp  and  current  components  can 
be  used  to  solve  for  the  horizontal  and  vertical 
components  of  points  along  horizontal  lines.  The 
restriction  of  remaining  somewhat  away  from  the 
center  and  rim  for  the  computed  results  will  still 
apply.  The  three  vertical  distances  from  the  ground 
plane  used  with  the  probe  are  now  used  when  solving 
for  theoretical  Jhorl2.  Jvert,  and  JtQtal.  The 

expression  for  total  current  would  still  be 

J total  *  <V  +  V>*  0 


where 


J  =  J  „„„  sin* 
p  p  max 


J,  =  J ,  cos* 

$  <p  max 

Now,  the  horizontal  and  vertical  components  are 
simply 

J,  =  J  sin*  +  J  cosd) 

horiz  p  Y  p 


y 


=  J  cos*  +  J  sin* 
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In  the  program  called  "DISK"  found  in 
Appendix  D,  it  was  convenient  to  transform  to 
rectangular  coordinates  and  iterate  along  the 
horizontal  limits  of  the  disk  at  a  particular 
height.  The  key  to  the  following  figures  is: 

J total  "  solid  llne  { - } 

Jhoriz  '  broken  line  ( - - *•) 

Jyert  -  dotted  line  (*••••) 


38  LAMBDA  DISK  MIDDLE  SCAN 


Figure  173:  Plot  of  current  components  for  the 
.38  wavelength  circular  disk  (uppe 
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Figure  174:  Plot  of  current  components  for  the  .50 
wavelength  circular  disk  (lower  scan) . 


th  circular  disk  (middle  scan) 


disk  (upper 


LAMBDA  DISK  MIDDLE  SCAN 


the  1.0  wavelength  circular  disk 
(upper  scan). 
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Relative  Current 


To  compare  the  validity  of  IR  detection  at 
this  stage  in  the  overall  project,  one  is  interested 
in  the  form  of  the  obtained  current  solution  rather 
than  a  closed  form  relationship  of  coating  current 
denisity  to  the  unperturbed  metal  plate  current 
density.  Thus,  for  the  .38  wavelength  square  plate 
and  disk,  a  comparison  follows  for  the  vertical 
magnetic  field  probe  results  which  correspond 
essentially  to  total  current  in  this  case,  the  IR 
results,  and  the  theoretical  solution  for  the  disk. 

The  technique  is  that  of  scaling  the  plots  from  this 
work  to  a  common  horizontal  distance  for  a  given 
horizontal  scan  and  scaling  the  vertical  distance 
from  the  zero  current  base  line  to  a  common  vertical 
distance  for  the  current  at  the  mid  point  of  the  scan. 
The  equipment  used  included  an  HP  9875A  Digitizer,  an 
HP  9835A  Minicomputer  and  an  HP  9872A  Plotter.  Thus, 
form  is  perserved  though  a  numerical  comparison  does 
not  apply. 

The  key  to  the  following  plots  is  as  follows: 
Magnetic  Field  Probe  -  smooth  line 
IR  Detection  Results  -  jagged  line 
Theoretical  Solution  -  xxxxx 


WAVELENGTH  CIRCULAR  DISK  (LOWER  SCAN) 
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lorrelation  of  theoretical,  probe,  and  IR  results  for 
he  .38  wavelength  circular  disk  (lower  scan). 


Horizontal  Position  (cm) 

Figure  184:  Correlation  of  theoretical,  probe,  and  IR  results  for 
the  .38  wavelength  circular  disk  (middle  scan). 
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Figure  185:  Correlation  of  theoretical,  probe,  and  IR  results  for 
the  .38  wavelength  circular  disk  (upper  scan). 
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Figure  188:  Correlation  of  the  magnetic  field  probe  and  IR  results 
for  the  -38  wavelength  square  plate  (upper  scan). 
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Now,  to  compare  the  square  plate  theoretical 
results  one  looks  at  the  bisecting  lines  of  the  1.0 
wavelength  square  plate  solved  by  Wilton  and  notes 
once  again  that  the  dominant  current  is  vertical  so 
that  the  vertical  magnetic  probe  measurement  applies 
though,  unfortunately,  is  the  least  symmetric  probe 
result  obtained.  The  1.0  wavelength  and  .50  wavelength 
square  plates  are  included  for  the  lower  scan  for 
plots  on  carbon/paraffin  over  metal  and  for  aquadac 
on  plexiglass.  Note  that  the  shapes  are  essentially 
the  same  but  the  variations  from  maximum  to  minimum 
current  density  levels  are  higher  for  aquadac  on 
plexiglass.  The  key  for  these  four  plots  remains 
as  in  the  previous  ones.  The  last  plot  is  a  comparison 
of  the  theoretical  solutions.  Wilton's  is  the  solid 
line  and  the  Body  of  Translational  results  from 
Appendix  C  are  "X"s.  Note  the  decrease  at  the 
midpoint  of  the  "X"  plot. 
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elation  of  theoretical,  probe,  and  IR  results  for  the 
wavelength  square  plate  (carbon/paraffin  on  metal) . 


0  WAVELENGTH  SQUARE  PLATE  (LOWER  SCAN) 
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elation  of  theoretical,  probe,  and  IR  results  for  the 
wavelength  square  plate  (carbon/paraffin  over  metal) . 


50  WAVELENGTH  SQUARE  PLATE  (LOWER  SCAN) 


re  192:  Correlation  of  theoretical,  probe,  and  IR  results  for 

the  .50  wavelength  square  plate  (aquadac  on  plexiglass) 


Horizontal  Position  (cm) 

Figure  193:  Theoretical  comparison  of  Body  of  Transluation 

Computer  Code  for  the  1.0  wavelength  square  plate. 


CHAPTER  VIII 


CONCLUSIONS 

The  qualitative  studies  on  the  infrared  detec¬ 
tion  of  microwave  induced  surface  currents  has  been 
extended  to  a  quantitative  investigation  for  flat 
plates.  Though  high  resolution,  infrared-obtained 
current  density  distributions  were  not  achievable 
with  the  present  state  of  coating  development  and 
available  equipment,  a  methodology  for  verifying 
the  IR  results  against  classic  approaches  was 
demonstrated. 

A  laboratory  facility  was  designed  and 
constructed  to  test  the  merits  of  the  IR  approach 
using  simple  shapes.  A  plane  wave  was  established  at 
the  target  and  verified  with  a  magnetic  field  probe. 

The  vertical  current  component  was  correlated  through 
magnetic  field  probe  measurements  with  the  theoretical 
prediction  for  the  perfectly  conducting  case.  The 
vertical  probe  measurements  for  dominant  vertical 
component  cases  detected  total  current  and  compared 
favorably  with  similar  patterns  seen  by  the  IR  on  both 
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carbon/paraffin  over  metal  and  aquadac  on  plexi¬ 
glass  . 


The  case  of  the  .38  wavelength  square  plate 
and  circular  disk  exhibited  a  high  degree  of 
correlation  with  theoretical  and  magnetic  field  probe 
results.  The  asymmetries  and  lack  of  high  resolution 
of  some  of  the  other  targets  could  be  attributed  in 
part  to  coating  inhomogeneities  and  separation  distance 
from  the  metallic  plate  as  well  as  a  digitization  and 
storing  process  that  eliminates  one  half  of  the 
available  data.  The  thermal  reservoir  effect  of  a 
large  aluminum  ground  plane  is  also  detrimental  for 
optimum  current  detection  based  on  a  heating  pattern. 

Supplementary  investigations  such  as  the 
angular  dependence  of  emissivity  and  subsequent 
utilization  of  a  correction  matrix  point  to 
applicability  for  complex  structui  where  probe 
measurements  become  increasing  more  difficult  and 
time  consuming. 

Having  developed  an  experimental  procedure 
to  relate  isotherm  unit  determinations  to  temperature, 
and  finally  current  density,  one  can  proceed  with 
further  study  of  the  coating  problem.  Verification 
should  continue  along  this  established  path  of  correla¬ 
tion  with  magnetic  field  probes  and  theoretical 


results  until  coatings  are  improved  and  simple 
shapes  have  been  studied.  As  target  complexity 
increases,  the  rapid  method  of  information 
collection  over  the  entire  object  using  infrared 
techniques  should  prove  to  be  increasingly 
significant  when  contrasted  with  classic  point 
by  point  measurements  and  mathematical  modeling 
techniques . 
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APPENDIX  A 

WAVE  DEMONSTRATOR 

A  ripple  tank  was  designed  and  constructed  for 
two-'dimensional  modeling  and  visualization  of  wave 
phenomena.  The  water  wave  interactions  should  roughly 
correspond  to  electromagnetic  wave  interactions  in  an 
anechoic  chamber.  This  appendix  develops  the  basic 
theory  of  water  wave  propagation  and  includes  data 
from  experiments  conducted  to  characterize  water  wave 
propagation  in  the  tank. 

Tank 

The  wave  demonstrator  tank  was  constructed  of 
a  5'  x  8*  sheet  of  3/16  inch  plexiglas,  with  2"  x  6" 
wood  sides  waterproofed  with  epoxy  paint  and  fastened 
with  1/4  inch  bolts  at  the  corners.  The  plexiglass 
sheet  was  fitted  into  a  3/4  inch  slot  around  the  wood 
sides  and  sealed  with  silicon  adhesive.  Figure  A1 
shows  the  construction  drawings.  The  assembled  wave 
generator  apparatus  is  seen  in  Figure  A2 .  The  center 
support  was  usually  necessary  to  prevent  sagging  of 
the  plexiglass  sheet  as  the  tank  could  be  filled  with 
up  to  3  cm  depth  of  water. 
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A  28  volt  D.C.  motor  turned  the  wave  generator 
drive  shaft,  which  rotated  two  eccentric  circular 
cams.  Ten  inch  aluminum  rods  set  on  the  cams,  which 
drove  the  probes  which  extended  into  the  water  and 
produced  the  water  wave .  The  circular  cams  were 
drilled  3/16  inch  off  center,  which  produced  a  near-? 
sinusoidal  motion  at  the  cam  of  3/32  inch  amplitude. 
The  pivot  ends  of  the  aluminum  rods  and  the  probe 
attachments  were  moved  toward  or  away  from  the  cams  to 
make  the  probe  oscillation  amplitude  variable  between 
3/32  inch  and  1/2  inch. 

Several  different  probes  were  attached  to  the 
wave  generator.  With  the  cams  rotated  in  phase,  a 
cylindrical  rod  was  attached  to  produce  plane  waves. 

A  single  round  bead  was  attached  to  one  rod  to  produce 
spherical  waves.  Two  element  antenna  array  simulation 
was  accomplished  with  a  two  bead/two  rod  arrangement. 
The  point  sources  could  also  be  driven  out  of  phase 
by  rotation  of  the  cams  on  the  drive  shaft. 

The  wave  generator  was  designed  to  be  clamped 
to  the  side  of  the  wave  demonstrator  tank,  or  it  may 
be  suspended  above  the  tank  for  experiments  away  from 
the  tank  sides. 

When  a  variable  voltage  power  supply  was  used, 
the  frequency  of  the  wave  generator  was  variable  up 
to  six  cycles  per  second.  The  maximum  voltage  which 
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could  be  applied  to  the  D.C.  motor  was  28  volts.  The 
graph  (Figure  A6)  shows  the  frequency  for  the  wave 
generator  as  a  function  of  applied  voltage,  and  was 
used  to  obtain  desired  frequencies  in  the  wave 
demonstrator  tank. 

Miscellaneous  Equipment 
A  point  light  source  suspended  from  the  ceil¬ 
ing  was  necessary  to  project  the  images  of  the  water 
waves  onto  the  floor  beneath  the  tank.  Also,  a  strobe 
light  was  used  to  stop  or  slow  the  apparent  motion  of 
the  v  xves  for  detailed  study  or  wavelength  measurement. 
A  white  screen  ruled  in  20  cm  x  20  cm  squares  was  used 
on  the  floor  beneath  the  tank  to  facilitate  viewing 
and  analysis  of  the  wa/es.  Black  and  white  photographs 
of  the  wave  images  on  the  floor  were  also  taken  by 
mounting  the  flash  above  the  tank  and  focusing  the 
camera  on  the  floor  from  the  side  of  the  tank.  The 
proper  camera  f-stop  for  tank  photographs  with  400 
ASA  film  was  f/8.5.  Absorbing  materiel,  including 
cloth  strips,  foam  rubber,  and  sloping  barrier 
inserts  could  be  added  to  the  tank  wherever  reflections 
were  not  desired. 

Theory  of  Water  Wave  Propagation 
The  following  analysis  is  from  Water  Waves  by 
N.  F.  Barber  (  1  ) .  A  more  detailed  derivation 
including  water  viscosity  and  surface  tension  may  be 
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found  in  Ripple  Tank  Studies  of  Wave  Motion  by  W. 
Llowarch  (2)  ,  and  a  circulation  approach  by  complex 
analysis  is  developed  in  Water  Waves  by  J.  J.  Stoker 
(3) •  The  approach  taken  in  this  appendix  is 
sufficient  to  understand  wave  propagation  for  ripple 
tank  modeling,  and  gives  results  which  correspond 
closely  to  those  observed  in  the  tank. 

The  motion  of  water  molecules  in  water  waves 
is  both  transverse  and  longitudinal — the  actual 
motion  is  approximately  elliptical.  The  molecules  at 
the  tops  of  the  wave  crests  move  forward,  in  the 
troughs  they  move  opposite  the  direction  of  the  wave, 
and  in  the  leading  and  trailing  slopes  the  molecules 
move  up  and  down,  respectively. 


Figure  A3.  Wave  Motion. 

The  motion  of  any  single  molecule  may  be 
broken  down  into  two  components:  Traverse  and 
longitudinal.  These  amplitudes,  A  and  B,  respectively, 
are  shown  in  Figure  A4 . 


Figure  A4 .  Wave  Amplitudes. 

Note  that  A  is  not  necessarily  equal  to  B,  and  in  a 
shallow  wave  tank#  A  will  generally  be  much  smaller 
than  B.  Also,  the  amplitudes  of  the  motion  A  and  B 
will  be  much  smaller  than  the  wave  length  A.  If  the 
shape  of  the  wave  is  assumed  to  be  sinusoidal,  then 
the  angle  between  the  vertical  and  the  normal  to 
the  water  surface  where  the  water  slope  is  the 
greatest  is  ,  shown  in  Figure  A5,  below. 


256 


The  ratio  of  the  horizontal  acceleration  of  the  water 
at  that  point  to  the  acceleration  of  gravity  is 
simply  the  tangent  of  =~^  ,  and  since  A<<X,  can  be 

2  it  A 

considered  equal  to  — .  Thus, 

Horizontal  Acceleration  (Max)  2ttA 
Gravitational  Acceleration  "T- 

and  thus 

Horizontal  Acceleration  (Max)  =  ?~~3.  . 

A 


Knowing  the  longitudinal  amplitude  of  the  sinusoid  B 
and  the  horizontal  acceleration,  the  angular  frequency 
oo  may  be  deduced: 

Horizontal  acceleration  (Max)  =  w2B, 
thus  =  (jJ2b 


so 


00 


/2vAq 
Y  BX 


2  TT 


The  period  T  equals  —  ,  so 


x  = 


gT2  a 


2  it  B 

Since  the  wave  velocity  V  equals  ~  , 

=  /nr 

Y  B  2  IT 


V 


(Al) 


Equation  Al  shows  how  the  wave  velocity  depends  on  the 
wavelength,  and  thus  the  frequency.  As  the  frequency 
increases,  the  wavelength  decreases,  and  thus  the  wave 
velocity  decreases.  i 

Also  as  the  water  gets  shallower,  the  ratio 
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—  in  equation  A1  gets  correspondingly  lower.  There- 
fore,  as  the  depth  of  the  water  increases  or 
decreases,  the  velocity  of  the  waves  in  turn  increases 
or  decreases. 

For  very  shallow  water  (for  example  5  mm  in 
the  wave  demonstrator  tank)  the  depth  of  the  water  D 
is  much  less  than  the  wavelength  X,  If  the  vertical 
velocity  component  of  the  water  is  assumed  to  also 

a 

follow  a  sinusoidal  function,  then  the  ration  —  equals 

2*D 
A  ’ 

Substituting  into  equation  A1 : 

V  =  v/~gD  (A2) 

Therefore,  for  very  shallow  water  (where  D<<A) ,  the 

wave  velocity  is  only  dependent  on  the  depth  of  the 

water,  and  little  dispersion  of  the  different  frequency 

components  oi  a  wave  should  occur. 

Comparing  the  predicted  results  from  equation 

A1  to  actual  results  in  the  wave  demonstrator  tank: 

A  m 

For  —  assumed  equal  to  1,  X  =  0.1m,  g  =  9.81  — 

33  s2 


The  experimental  value  tank  measurements  was 

V  =  0.5  —  .  For  a  water  depth  of  5  mm,  using  equation 
s 


A2: 


V  =  >/gD* 

=  0.7  “ 
s 

The  experimental  value  from  tank  measurements  was 

V  =  0.4  —  .  Therefore,  the  theoretical  results 
s 

corresponded  favorably  to  the  actual  results  obtained 
from  the  wave  demonstrator  tank. 

Theory  of  Wave  Demonstrator  Operation 
The  wave  demonstrator  tank  provided  a  visual¬ 
ization  in  two  dimensions  of  wave  phenomena  by 
comparison  with  water  waves.  The  water  waves  exhibit 
many  of  the  characteristics  of  electromagnetic  waves: 
They  reflect  at  their  incident  angle  with  a  180°  phase 
change;  they  interfere  constructively  and  destructively 
to  produce  standing  waves  and  interference  patterns; 
they  travel  slower  through  shallower  water  to 
simulate  higher  indices  of  refraction,  etc. 

Although  the  waves  themselves  were  very 
difficult  to  see  in  the  tank,  they  were  easily  observed 
by  projecting  their  images  onto  the  floor  beneath  the 
tank.  When  a  "point"  light  source  was  suspended  above 
the  tank  in  a  dark  room  and  a  sheet  of  white  paper  was 
placed  beneath  the  tank,  the  convex  and  concave 
curvature  of  the  waves  acted  like  converging  and 
diverging  lenses,  projecting  bright  and  dark  lines  on 
the  paper.  The  arrangement  is  shown  in  Figure  A2, 


which  also  shows  the  distances  from  the  tank  to  the 
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light  and  the  paper  used  for  the  experiments  described 
in  this  appendix.  These  distances  were  not  optimum 
for  all  wavelengths  and  water  depths,  but  were  con¬ 
venient  and  suitable  for  a  general  range  of  use. 

Experimental  Wave  Characterization 
Plane  Waves 

The  wave  generator  was  clamped  to  the  tank 
side  and  the  plane  wave  generator  attachment  (a  4' 
wooden  dowel)  was  connected  to  the  wave  generator  arms. 
Using  a  raeterstick  and  a  stopwatch,  measurements  of 
wave  velocity  and  wavelength  were  made  at  various 
water  depths  from  5  mm  to  30  mm  and  frequencies  rang¬ 
ing  between  1.5  and  5.5  cycles  per  second.  The 
experimental  data  and  calculations  of  wave  velocity, 
frequency,  wavelength,  and  voltage  applied  to  the 
motor  are  available  in  a  laboratory  notebook  at  the 
Air  Force  Academy.  The  data  is  shown  in  Figures  A6 
through  A8  as  graphs  of  frequency  as  a  function  of 
voltage  applied  to  the  motor,  and  wave  velocity  and 
wavelength  of  functions  of  frequency  for  various 
water  depths.  These  graphs  were  used  to  choose  the 
required  water  depth  and  motor  voltage  to  produce  a 
desired  frequency,  wavelength,  and/or  wave  velocity. 
Circular  Waves 

Because  the  circular  wave  probe  (a  wooden 
bead)  displaces  much  less  water  than  the  plane  wave 
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Figure  A7:  Velocity  vs  Frequency 
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probe,  and  due  to  the  j  dependence  of  the  amplitude 

of  a  two  dimensional  wave  and  the  great  attenuation 
over  distance  of  a  water  wave,  the  images  formed  of 
a  circular  wave  were  considerable  dimmer,  especially 
at  low  frequency  and  shallow  water  depth. 

Consequently,  it  was  recommended  that  a  minimum  water 
depth  of  10  mm  or  minimum  wave  generator  frequency  of 
3  cycles  per  second  be  used  to  obtain  easily  visible 
circular  wave  images.  The  included  calibration 
charts  were  used  to  set  approximate  motor  voltages  and 
water  depths  for  desired  frequencies,  wavelengths,  and 
wave  velocities. 

Visual  Characteristics 

The  waves  in  the  tank  projected  sharper  images 
on  the  floor  as  the  water  depth  was  increased.  How¬ 
ever,  dispersion  of  the  higher  frequency  components 
becomes  more  of  a  problem  with  increased  water  depth, 
as  shown  in  the  wave  velocity  vs.  frequency  chart  in 
Figure  A  7 .  The  best  balance  between  image  clarity 
and  reduced  dispersion  was  determined  experimentally 
to  occur  at  approximately  15  mm  water  depth,  for  wave¬ 
lengths  associated  with  antenna  modeling  and  acoustic 
characterization  experiments.  This  is  especially 
important  for  circular  waves  produced  by  a  spherical 
probe,  since  the  amplitude  of  the  plane  waves  produced 
by  a  cylindrical  probe  is  greater.  For  reflection. 
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refraction,  and  interference  experiments  where  higher 
frequencies  (5  to  6  cycles  per  second)  were  used, 
slower  waves  with  less  dispersion  were  obtained  and 
seen  easily  using  shallower  water  depths  (less  than 
10  mi.  . 

Operation  of  the  Wave  Demonstrator 
Water  depth  and  wave  generator  frequency  most 
suitable  for  proper  wavelength  and  velocity  and 
sharpest  viewing  were  selected.  For  short  wavelengths 
with  maximum  clarity,  a  shallow  water  depth  (10  mm  or 
less)  is  most  desirable,  providing  the  least  dispersion 
and  low  wave  velocity  at  the  shortest  wavelengths. 

For  demonstrations  of  varying  wave  velocity  with  vary¬ 
ing  frequency  (dispersion),  a  deeper  water  depth  (30 
mm)  were  most  effective.  Where  waves  were  scaled  to 
the  dimensions  of  an  anechoic  chamber  containing  a 
reflector  and  scattering  object  a  compromise  between 
wave  image  brightness  and  frequency  dispersion  was 
achieved  with  a  water  depth  of  about  15  mm.  Once  the 
proper  water  depth  had  been  selected,  the  desired 
wavelength  or  wave  velocity  was  found  from  the  included 
graphs.  After  selecting  the  required  wavelength  or 
wave  velocity  from  the  family  of  water  depth  curves, 
the  required  frequency  was  read  from  Figure  a8 •  The 
applied  voltage  for  the  motor  was  determined  from  the 
frequency  vs.  voltage  graph. 


The  wavelengths  in  the  graphs  were  measured 
from  the  images  projected  on  the  paper  by  the  setup 
in  Figure  A2 .  To  convert  these  lengths  to  actual 
wavelengths  of  the  water  waves  in  the  tank,  the 
dimensions  were  multiplied  by  — the  ratio  of 

distances  of  the  light  from  the  tank  and  floor. 

With  the  tank  filled  to  the  required  depth, 
the  light  was  suspended  approximately  1.7  meters  above 
the  water,  and  a  white  screen  positioned  approximately 
one  meter  below  the  tank.  For  the  best  wave  images, 
the  height  of  the  probe  was  adjusted  so  that  the 
probe  comes  as  close  as  possible  to  the  plexiglass 
bottom  of  the  tank,  without  touching  the  bottom  at 
the  lowest  point  in  its  motion,  or  breaking  contact 
with  the  water  surface  at  the  top.  The  best  images 
usually  occurred  with  minimum  wave  generator  amplitude, 
except  when  using  the  circular  wave  generator  at  low 
frequency,  where  a  larger  amplitude  will  be  needed. 

The  images  of  the  waves  in  the  tank  were  clearly 
visible  on  the  screen  beneath  the  tank  when 
illuminated  by  a  "point"  light  source  above  the  tank. 

The  antenna  patterns  from  corner  reflector 
antennas  were  modeled  to  scale  in  the  tank  using  the 
circular  wave  generator  probe  and  placing  straight 
rigid  barriers  at  an  angle  behind  the  probe. 
Recommended  dimensions  for  corner  reflectors  were 


found  in  "Gains  of  Finite-Size  Corner  Reflector 
Antennas"  by  H.  V.  Cottony  and  A-  C.  Wilson  (4), 
Parabolic  transmitting  antennas  were  modeled  in  the 
tank  by  generating  circular  waves  at  the  focus  of  the 
parabolic  reflector.  A  waveguide  and  microwave  horn 
antenna  were  also  modeled  in  the  tank  by  propagating 
plane  waves  down  a  ^  wide  constriction  in  the  tank 

and  into  a  diverging  set  of  rigid  barriers.  Recom¬ 
mended  dimensions  for  horn  antennas  may  be  found  in 
Antennas  by  John  D.  Kraus  (5). 

The  photographs  that  follow  were  taken  using 
only  a  flash  mounted  above  the  tank  in  an  otherwise 
dark  room.  An  open  shutter  at  f/8.5  was  used  with 
400  ASA  film.  Note  the  relatively  plane  wavefronts 
after  approximately  three  divisions  on  the  screen. 

The  dark  vertical  lines  are  meter  stick  barriers 
modeling  the  sides  of  "perfectly"  reflecting  anechoic 
chamber  walls.  Corner  reflector  aperature  angles  and 
antenna  to  apex  distances  were  varied  to  investigate 
gain  and  plane  wavefront  parameters. 


Figure  A10 :  9 
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APPENDIX  C 

BODY  OF  TRANSLATION  SQUARE  PLATE  RESULTS 

The  square  plate  problems  were  run  with  a 

Bodies  of  Translation  Code  (BOT)  called  A-STAR.  This 
code,  developed  by  McDonnell  Douglas  Research  Labs, 
models  bodies  which  have  a  constant  cross  sectional 
geometry  such  as  flat  plates  and  cylinders.  The  code 
predicts  two  orthogonal  components  of  surface  current 
distribution  and  radiated  near  and  far  electric  and 
magnetic  fields.  The  code  also  has  the  capability  of 
modeling  enclosing  encaps  on  the  geometry  as  well  as 
attached  thin  wires. 

The  concept  of  this  model  assumes  model  current 
distributions  of  the  form  exp(jntA/L)  along  the  length 
of  the  body.  This  model  is  presently  limited  to 
bodies  of  about  10  square  wavelengths. 
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.38  LAMBDA  SQUARE  PLATE 

Figure  C6:  Normalized  current  contour  plot  for  the 
.38  wavelength  square  plate  (BOT) . 
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.50  LAMBDA  SQUARE  PLATE 


Figure  C7 :  Normalized  current  contour  plot  for  the 

.50  wavelength  square  plate  (BOT). 


1.3  LflMBDR  SQURRE  PLATE 


.a 


Normalized  current  contour  plot  for  the 
1.3  wavelength  square  plate  (BOT) . 


Table  Cl 


Current  Components  from  the  Body 
of  Translation  Code  for  a  1.0 
wavelength  square  plate 


^horizontal  (A/m) 
. 763E-09 
. 226E-09 
, 864E-09 
. 355E-09 
. 115E-08 
. 502E-09 
. 135E-08 
. 455E-11 


661E-09 


472E-02 
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APPENDIX  D 

COMPUTER  PROGRAMS 

CAMCAL  -  Calibration  for  Thermovision 
J  vs  T  -  Quadratic  fit  from  empirical  data 
IRSTAT  -  Statistical  Analysis  on  the 

Initial  and  Final  IR  data  then 
stores  it 

MEGAIR  -  Calculates  temperature  and 

current  from  data  -  also  graphs 
PDISC  -  Bouwkamp  disk  solution 
DISK  -  Cartesian  and  Total  current 


calculation 
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zontal  Magnetic  Field  Frobe  Measurements  for  the 
8  wavelength  circular  disk. 


.38  wavelength  square  plate. 
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Figure  E7:  Horizontal  Magnetic  Field  Probe  Measurements  for  the 

1,3  wavelength  circular  disk. 


‘Frequency  2.45  GHz 


Horizontal  Magnetic  Field  Probe  Measurements  for  the 
1.3  wavelength  square  plate. 
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APPENDIX  F 

EMISSIVITY  DEPENDENCE  ON  ANGLE 

Emissivity  is  the  absorption/emission 
properties  of  a  substance  as  compared  to  the  isotropic 
and  perfectly  absorbing/emitting  properties  of  a 
black  body  radiator.  To  determine  the  true  temperature 
of  an  elemental  area  on  an  object  surface  using  an 
infrared  system,  the  directional  nature  of  the  IR 
photon  emission  was  investigated.  A  simple  point  by 
point  correction  using  a  stored  object  matrix  taken 
in  an  ambient  temperature,  non-illuminated  state  was 
found  to  be  effective  with  temperature  elevations  up 
to  5°  K  above  ambient.  This  work  was  reported  through 
presentations  and  publications  to  SPIE  -  The  Inter¬ 
national  Society  for  Optical  Engineering  (  4  )  and 
the  IEEE  (  5 ) •  It  represents  an  important  gain  in 
knowledge  that  is  most  applicable  for  viewing  angles 
that  are  other  than  normal  incidence. 

Experimentally  Determined  Emissivities 
In  determining  the  surface  temperature  distribu¬ 
tion  of  three  dimensional  objects  by  thermographic 
analysis,  the  problem  of  directional  emissivity — 


angular  and  material  dependence  of  infrared  emission 
from  varying  surfaces — must  be  accounted  for.  At  a 
given  temperature,  different  materials  emit  differing 
intensities  of  infrared  radiation,  and  the  intensity 
is  also  directionally  dependent.  This  could  cause 
problems  with  correlating  the  thermographic  output 
from  an  infrared  camera  with  actual  surface  tempera¬ 
tures  on  objects  of  various  material  coatings  or 
complex  geometries.  A  simple  experiment  was  conducted 
using  a  hollow  aluminum  sphere  coated  with  a  carbon 
base  paint  and  maintained  at  a  constant  temperature 
with  an  interior  water  bath  held  at  the  boiling  point. 
The  arrangement  is  seen  in  Figure  (Fl )  .  The  Thermo¬ 
vision,  depicted  the  sphere  in  ten  colors,  each 
representing  1°  C,  as  nine  concentric  rings  around  a 
center  circle,  Figure  (f2) •  The  Thermovision  depiction 
of  a  ten  degree  temperature  range  contrasted  with  the 
nearly  uniform  temperature  profile  actually  present. 
This  demonstration  clearly  indicated  the  requirement 
to  correct  for  geometry  to  obtain  accurate  results. 

The  influence  of  the  emissivity  dependence  on  thermo- 
graphically  determined  temperature  distributions  was 
also  noted  by  Hsieh  and  Ellington  (  2  )  in  1977. 
Idealized  plots  of  emissivity  and  angular  dependence 
can  be  found  in  several  sources  (6  ).  However,  because 
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various  coating  preparations  were  to  be  used  in  the 
thermographic  analysis,  it  was  felt  necessary  to 
experimentally  determine  the  directional  emissivity 
of  various  potential  coating  materials. 

For  this  determination,  an  apparatus  to 
determine  emissivity  at  various  viewing  angles  for 
several  different  materials  and  coatings  was  built, 
and  is  shown  in  Figure  (F3).  Patches  of  different 
resistive  coatings  which  could  be  used  in  our 
analysis  were  applied  to  the  side  of  a  water-filled 
aluminum  box.  A  15  cm  diameter  spherical  aluminum 
chamber  coated  inside  with  lampblack  and  welded  into 
the  side  of  the  box  with  a  5  cm  diameter  circular 
opening  at  the  front  of  the  box  served  as  a  simulated 
blackbody  radiator,  with  which  the  emissivities  of 
the  various  coating  patches  could  be  compared.  The 
emissivity  of  the  blackbody  reference  was  taken  as 
er=l,  and  verified  to  be  isotropic  by  viewing  the 
blackbody  with  the  Thermovision  camera  at  angles 
between  0°  and  70°  to  the  normal  of  the  face  of  the 
box  (the  detected  IR  radiation  from  the  blackbody 
simulator  did  not  vary  with  angle  until  the  camera 
could  not  longer  see  inside  the  lip  of  the  simulator 
opening  at  70°  from  the  normal) .  The  water  was 
maintained  at  boiling  (about  92°  C) ,  and  thermographic 


(a)  Experimenta 
Lgure  F3:  Black 
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data  was  taken  at  angles  between  0°  and  approximately 
90°  from  the  normal  to  the  box  face  with  an  AGA 
Thermovision  680  interfaced  with  a  HP  9845B  mini¬ 
computer  described  in  the  previous  section.  The 
emissivities  as  a  function  of  angle  were  calculated 
using  the  following  formula,  explained  in  the  AGA 
Thermovision  680  Operating  Manual  ( 1  ) : 


e(e) 


4i+er(Ir-Ia» 


a 


where  Ai  is  the  difference  in  isotherm  levels  on  the 
Ihermovision  between  the  test  material  (varying  with 
angle  0)  and  the  blackbody  reference,  and  Ir,  I  ,  and 
IQ  are  the  isotherm  levels  of  the  blackbody  reference, 
ambient,  and  test  material  temperatures,  respectively. 

Figure  (F4)  shows  the  position  of  the  material 
on  the  simulator  face  and  a  digital  printout  at  10° 
from  normal  incidence.  Only  the  center  portion  was 
used  to  determine  coverage  isotherm  readings  used  in 
the  calculations. 


Figure  (f5)  depicts  experimental  results  from 


the  emissivity  test  apparatus  for  teledeltos  paper  and 
smooth/rough  surface  preparations  of  aquadac. 
Teledeltos  is  a  resistive  paper,  and  aquadac  is  also 
available  resistive  coating,  discussed  earlier.  Also 
shown  for  comparison  are  the  results  for  a  typical 
dielectric  (wood)  and  aluminum,  on  which  the  samples 
were  attached.  Curves  for  wood  and  aluminum  are 
taken  from  the  literature.  The  dielectric  nature  of 
the  coating  appeared  to  account  for  the  decreasing  IR 
emittance  from  the  normal  to  tangential  viewing  angles 


Figure  F5:  Total  Directional 

Emittance  Distributions 


It  was  interesting  that  roughness  appeared  to 
be  a  factor  for  material  emissivity  and  this  prompted 
a  brief  Scanning  Electron  Microcrope  observation  of 
the  coating  surfaces  to  be  used  in  the  IR  detection 
of  surface  currents.  The  results  are  seen  in  Figures 
(F6 )  and  (F7). 


(a)  Carbon/paraffin  (1300X) 


Figure  F7: 


(b)  Carbon/paraffin  (5200X) . 

Micrographs  of  the  carbon/paraffin 
coating. 


Correction  for  Angular 
Dependency  of  Emissivrty 
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As  was  discussed,  a  correction  was  necessary 
since  an  accurate  determination  of  surface  temperature 
was  a  prerequisite  for  obtaining  surface  current 
information.  The  emissivity  variation  with  viewing 
angle  was  determined  to  be  the  needed  correction.  The 
theoretical/  mathematical  approach  was  presented  (3  ) 
and  later  replaced  by  a  simpler  experimental  one, 
summarized  here.  For  temperature  variations  that 
typically  do  not  exceed  5°  C  it  was  found  that  a 
matrix  of  temperature  values  taken  of  the  object  at 
a  steady  state  ambient  temperature  condition  prior  to 
irradiation  by  microwave,  could  be  used  effectively 
for  emissivity  correction  for  geometry.  Experimental 
verification  of  the  validity  of  the  correction  for 
angular  dependence  of  emissivity  by  simple  subtraction 
of  the  correction  matrix  from  a  final  "heated" 
temperature  matrix  was  carried  out  over  a  range  of 
temperatures  above  ambient  for  the  hollow  sphere. 

The  procedure  for  determining  currents  on 
complex  shapes  would  involve  storage  of  a  correction 
matrix  obtained  from  an  average  of  several  frames 
taken  of  the  object  at  ambient  temperature.  A 
several  frame  average  of  an  illuminated  or  "final" 
matrix  is  stored.  The  correction  matrix  is  then 


subtracted  from  the  illuminated  to  yield  a  difference 
matrix  in  ISU  units  which  now  has  accounted  for 
emissivity  variation  with  angle.  The  standard 
conversions  to  temperature  and  current  are  then 
made.  The  following  five  figures  depicts  this 
procedure  applied  to  an  early  aquadac/plexiglass 
flat  plate. 


Horizontal  <X)  Axis 


Figure  F10:  Difference  Matrix  (ISU) 


Horizontal  OO  Axis 
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